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Abstract 

A partition (Ci, C 2 , Cq) of G = {V, E) into clusters of strong (respectively, weak) diameter d, such 
that the supergraph obtained by contracting each Ci is Gcolorable is called a strong (resp., weak) (d, t)- 
network-decomposition. Network-decompositions were introduced in a seminal paper by Awerbuch, 
Goldberg, Luby and Plotkin in 1989. Awerbuch et al. showed that strong {exp{0{i/\og n log log n)}, 
ea;p{0(-\/log n log log n)})-network-decompositions can be computed in distributed deterministic time 
exp{0 {\/log n log log n)}. Even more importantly, they demonstrated that network-decompositions 
can be used for a great variety of applications in the message-passing model of distributed computing. 

The result of Awerbuch et al. was improved by Panconesi and Srinivasan in 1992: in the latter re¬ 
sult d = £ = exp{0{y/log n)}, and the running time is exp{0{i/log n)} as well. In another remarkable 
breakthrough Linial and Saks (in 1992) showed that weak (O(logn), 0(log n))-network-decompositions 
can be computed in distributed randomized time 0(log^ n). Much more recently Barenboim (2012) de¬ 
vised a distributed randomized constant-time algorithm for computing strong network decompositions 
with d = 0(1). However, the parameter £ in his result is 0(n^/^+*^). 

In this paper we drastically improve the result of Barenboim and devise a distributed randomized 
constant-time algorithm for computing strong (0(1), 0(n'^))-network-decompositions. As a corollary 
we derive a constant-time randomized 0(n'^)-approximation algorithm for the distributed minimum 
coloring problem, improving the previously best-known 0(n^/^“'"‘^) approximation guarantee. We also 
derive other improved distributed algorithms for a variety of problems. 

Most notably, for the extremely well-studied distributed minimum dominating set problem cur¬ 
rently there is no known deterministic polylogarithmic-time algorithm. We devise a deterministic 
polylogarithmic-time approximation algorithm for this problem, addressing an open problem of Lenzen 
and Wattenhofer (2010). 


1 Introduction 

1.1 Network-Decompositions 

In the distributed message-passing model a communication network is represented by an n-vertex graph 
G = {V, E). The vertices of the graph host processors that communicate over the edges. Each vertex has 
a unique identity number (ID) from the range {1,2, ...,n}. We consider a synchronous setting in which 
computation proceeds in rounds, and each message sent over an edge arrives by the beginning of the next 
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round. The running time of an algorithm is the number of rounds from the beginning until all vertices 
terminate. Local computation is free. 

A strong (respectively, weak) diameter of a cluster C C 1/ is the maximum distance distc(f;) (n, v) 
(resp., distG(u,u)) between a pair of vertices u,v G C, measured in the induced subgraph G{C) of C 
(resp., in G). A partition (Ci, 6 * 2 ,..., Cg) of G = {y,E) into clusters of strong (resp., weak) diameter 
d, such that the supergraph Q = (V,T), V = {Ci, (72,..., Cg}, £ = {{Ci,Cj) \ Gi,Cj GV,i G 

Ci,Vj G Gj,{vi,Vj) G E} obtained by contracting each Ci is Acolorable is called a strong (resp., weak) 
{d, i)-network-decomposition. 

Network-decompositions were introduced in a seminal paper by Awerbuch et al. [3]. The au¬ 
thors of this paper showed that strong {exp{O (Vlog ^ log log n)}, exp{0{-\/\og n log log n)})-network- 
decompositions can be computed in distributed deterministic exp{0{y/\og\og\ogn)} time. Even more 
importantly they demonstrated that many pivotal problems in the distributed message passing model 
can be efficiently solved if one can efficiently compute (d, t')-network-decompositions with sufficiently 
small parameters. In particular, this is the case for Maximal Independent Set, Maximal Matching, and 
(A -|- l)-Vertex-Coloring. 

The result of [3] was improved a few years later by Panconesi and Srinivasan |45j who devised a de¬ 
terministic algorithm for computing strong (exp{0{^ylogn)}, exp{0{y/\ogn) })-network-decompositions 
in exp{0{yJ\og n)} time. Awerbuch et al. [T] devised a deterministic algorithm for computing strong 
((7(log n), 0(log n))-network-decomposition in time exp{0{y/logn)}. Around the same time Linial and 
Saks [39] devised a randomized algorithm for computing weak (0(log n), 0(log n))-network-decompositions 
in 0{log^ n) time. More generally, the algorithm of |39] can compute weak (A, (9(n^/^ logn))-network- 
decompositions or weak ((9(n^/^), A)-network-decompositions in time 0(A • n^/^\ogn). 

Observe, however, that all these algorithms [3ll45[[39] require super-logarithmic time, for all choices of 
parameters. In ICALP’12 the first-named author of the current paper [5] devised a randomized algorithm 
for computing strong ((9(1), n^/^'’'^)-network-decomposition in (9(l/e) time. Unlike the algorithms of 
|3lll5l[39], the algorithm of |5] requires constant time. Its drawback however is its very high parameter 
i = In the current paper we alleviate this drawback, and devise a randomized algorithm for 

computing strong (exp{(9(A)}, n^/'^)-network-decomposition in time exp{0{X)}. In other words, the 
parameter A of our new decompositions can be made for an arbitrarily small constant e > 0, while 
the running time is still constant (specihcally, exp{0{l/e)'\). 

1.2 Constant-Time Distribnted Algorithms 

In their seminal paper titled ’’What can be computed locally?” [l3] Naor and Stockmeyer posed the 
following question: which distributed tasks can be solved in constant time? This question is appealing 
both from theoretical and practical perspectives. From the latter viewpoint it is justihed by the emergence 
of huge networks. The number of vertices in the latter networks may be so large that even mildest 
dependence of the running time on n may make the algorithm prohibitively slow. 

Naor and Stockmeyer themselves |l3| showed that certain types of weak colorings can be computed in 
constant time. A major breakthrough in the study of distributed constant time algorithms was achieved 
though a decade after the paper of [l3] by Kuhn and Wattenhofer [M]. Specifically, Kuhn and Wattenhofer 
|34j showed that an 0{-\/k/S}/'^ log A)-approximate minimum dominating se10 can be computed in 0{k) 
randomized time. Here A = A(G) is the maximum degree of the input graph G, and k is a positive 
possibly constant parameter. 

An approximation algorithm for another fundamental optimization problem, specifically, for the min¬ 
imum coloring problem, was devised by Barenboim [5] as an application of his aforementioned algorithm 
for computing network-decompositions. Specihcally, it is shown in |5] that an 0(n^/^^'^)-approximation 

^ A subset U C in a graph G = (U, E) is a dominating set if for every v £ V\U there exists u £ U, such that [u, v) £ E. 
In the minimum dominating set (henceforth, MDS) problem the goal is to find a minimum-cardinality dominating set of G. 
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for the minimum coloring problem can be computed in 0(l/e) randomized time. (In the minimum color¬ 
ing problem one wishes to color the vertices of the graph properly with as few colors as possible.) Observe 
that since approximating the minimum coloring problem up to a factor of is NP-hard |27l |2^ 149] . 
the algorithm of [5] inevitably has to employ very heavy local computations. 

In the current paper we employ our improved network-decomposition procedure to come up with a 
significantly improved constant-time approximation algorithm for the minimum coloring problem. Specif¬ 
ically, our randomized algorithm provides an 0(n^)-approximation for the minimum coloring problem 
in exp{0{l/e)'\ time, for an arbitrarily small constant e > 0. We also devise a randomized 0{rf)- 
approximation algorithm for the minimum t-spanner problem with running time exp{0{l/e)'\ -|- 0(f), 
for any arbitarily small constant e > 0. (A subgraph G' = {V,H) of a graph G = {V,E), H C E, is a 
t-spanner of G if for every u,v ^ V, distG'(^) v) < t ■ distG(?.i, v). In the minimum t-spanner problem the 
objective is to compute a t-spanner of the input graph G with as few edges as possible.) 

Ajtai et al. [2] showed that triangle-free n-vertex graphs admit an 0{^/nj \/logn)-coloring. This exis¬ 
tential bound was shown to be tight by Kim [29]. We devise a randomized 0(n^/^'’"'^)-coloring algorithm 
for triangle-free graphs with running time 0(l/e). More generally, we devise a randomized 
coloring algorithm for graphs of girth greater than g = 2k, k > 2, with running time 0(l/e^). Both results 
apply for any arbitrarily small e > 0, and, in particular, they show that such graph can be colored with a 
reasonably small number of colors in constant time. Together with our drastically improved constant-time 
approximation algorithm for the minimum coloring problem, these results significantly expand the set of 
distributed problems solvable in constant time. 

Most our algorithms for constructing network-decompositions use only short message^ (he., messages 
of size O(logn) bits), and employ only polynomially-bounded local computations. Although in general 
graphs our algorithms for 0(n^/’^)-approximate minimum coloring require large messages, our 0(re^/^"’“'^)- 
coloring and 0(n^/^"*“’^)-coloring algorithms for triangle-free graphs and graphs of large girth employ short 
messages. Hence the latter coloring algorithms are suitable to serve as building blocks for various tasks. 
Despite that the number of colors is superconstant, in many tasks it does not affect the overall running 
time, so the entire task can be performed very quickly. For example, if the colors are used for frequency 
assignment or code assignment tasks, the running time will not be affected by the number of colors. 
Instead, the range of frequencies or codes will be affected. However, this is unavoidable in the worst case, 
in view of the lower bounds on the chromatic number of triangle free graphs and graph of large girth. 
1.3 The Minimum Dominating Set Problem 

The MDS problem is one of the most fundamental classical problems of distributed graph algorithms. Jia 
et al. [28| devised the first efficient randomized 0(log A)-approximation algorithm for the MDS problem 
with running time 0(lognlog A). Their result was improved and generalized by Kuhn and Wattenhofer 
[34] who devised an 0(fc)-time randomized 0(\/fcA^/'^log A)-approximation algorithm for the problem. 

The results of |28l [34] spectacularly advanced our understanding of the distributed complexity of the 
MDS problem. However, both these algorithms [281 [34] are randomized, and no efficient deterministic 
distributed algorithms with a non-trivial approximation guarantee for general graphs are currently known. 
Lenzen and Wattenhofer m devised such algorithms for graphs with bounded arboricity. Below we 
provide a quote from their paper: 

”To the best of our knowledge, the deterministic distributed complexity of MDS approximation on general 
graphs is more or less a blind spot, as so far neither fast (polylogarithmic time) algorithms nor stronger 
lower bounds are known”. 

In this paper we address this blind spot and devise a deterministic 0(n^/^)-approximation algorithm for 
the MDS problem with time 0((logn)^“^). Similarly to our approximation algorithms for the minimum 
coloring and the minimum t-spanner problems, this algorithm is also a consequence of our algorithms 

^The only exceptions are weak network-decompositions from Section [6.31 
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for constructing network-decompositions. However, for the MDS we use a deterministic version of these 
algorithms, while for the minimum coloring and minimum t-spanner problems we use a randomized 
version. Also, we present a variant of our MDS approximation algorithm that employs only polynomially- 
bounded local computations, requires 0((logn)^“^) time, and provides an log A) approximation. 

1.4 Additional Results 

We also use our algorithms for computing network-decompositions for devising algorithms for computing 
low-intersecting partitions. Low-intersecting partitions were introduced by Busch et ah m in a paper 
on universal Steiner trees. A low-intersecting {a, (3,^)-partition P of a graph G is the partition of the 
vertex set V such that: (1) Every cluster C in V has strong diameter at most a ■ 7 . 

(2) For every vertex u G E, a ball B-^{y) of radius 7 around v intersects at most f3 clusters of V. 

Busch et al. showed that given a hierarchy of low-intersecting partitions with certain properties (see 
m for details) one can construct a universal Steiner tree. (See [Mj for the definition of universal Steiner 
tree.) Also, vice versa, given universal Steiner tree they showed that one can construct a low-intersecting 
partition. They constructed a low-intersecting partition with a = 4^,/? = k ■ and arbitrary 7 . 

We devise a distributed randomized algorithm that constructs low-intersecting (( 0 ( 7 )^, 7 )- 
partitions in time ( 0 ( 7 ))*^ log^'^^ n in general graphs and in ( 0 ( 7 ))^ • exp{0{-\/\og log n)} time in graphs of 
girth <7 > 6. This algorithm employs only short messages and polynomially-bounded local computations. 

Comparing this result with the algorithm of Busch et al. [13] we note that the partition of [T3| has 
smaller radius. (It is 7 • (0(1))^ instead of ( 0 ( 7 ))^ in our case.) On the other hand, the intersection 
parameter (3 of our partitions is smaller. (It is instead of k ■ n^/^.) In particular, the intersection 
parameter in the construction of [H] is always D(logn), while ours can be as small as one wishes. Finally, 
and perhaps most importantly, the algorithm of [TTj is not distributed, and seems inherently sequential. 

1.5 Comparison of Our and Previous Techniques 

Basically, our algorithms for computing network-decompositions can be viewed as a randomized variant 
of the deterministic algorithm of Awerbuch et al. [3|. The algorithm of Awerbuch et al. [3] computes 
iteratively ruling sets for subsets of high-degree vertices in a number of supergraphs. These supergraphs 
are induced by certain graph partitions which are computed during the algorithm. (A subset U O V 
of vertices is called an (a, /I)-ruling set if any two distinct vertices u,u' ^ U are at distance at least a 
one from another, and every v £ V \ U not in a ruling set has a ’’ruler” u £ U at distance at most 
f3 from V.) As a result of the computation the algorithm of |3| constructs a partition into clusters of 
diameter at most a, such that the supergraph induced by this partition has arboricity at most j3. The 
algorithm of [3| then colors this partition with 0{/3) colors in time 0(/31ogn) • 0(a). (The running time 
of the algorithm is 0(/31ogn) when running on an ordinary graph. The running time is multiplied by a 
factor of 0(a), because the coloring algorithm is simulated on a supergraph whose vertices are clusters 
of diameter 0(a).) The fact that the running time in the result of |3] is (roughly speaking) the product 
a • /? of the parameters of the resulting network-decomposition is the reason that Awerbuch et al [3] made 
an effort to balance these parameters, and set both of them to be equal to exp{0{ n log log n)}. 
The algorithm of Panconesi and Srinivasan [35] is closely related to that of [3] except that it invokes 
a sophisticated doubly-recursive scheme for computing ruling sets via network-decompositions, and vice 
versa. This ingenious idea enables [33] to balance the parameters and running time better. Specifically, 
they are all equal to 

Our algorithm is different from [1135] in two respects. First, we replace a quite slow (it requires 
O(logn) time) deterministic procedure for computing ruling sets by a constant-time randomized one. 
Note that generally computing (0(1), 0(l))-ruling sets requires D(log* n) time [38], but we only need to 
compute them for high-degree vertices of certain supergraphs. This can be easily done in randomized 
constant time. Second, instead of coloring the resulting partition with 0{/3) colors in 0(/31ogn) • 0(a) 
time, we color it in 0{f3 ■ n^) colors in 0(l/e) • 0(a) time by a simple randomized procedure, or in 
0(/3^1og*'*^ n) colors in 0{t) ■ 0(a) time, for a parameter t > 0, by a deterministic algorithm Arb-Linial 
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Hence the number of colors is somewhat greater than in OHS], but the running time is constant. 
The algorithm of Linial and Saks [39| is inherently different from both mm and from our algorithm. 
It runs for O(logn) phases, each of which constructs a collection of clusters of diameter O(logn) at 
pairwise distance at least 2 which covers at least half of all remaining vertices. The running time of the 
algorithm of [39], similarly to [3] and |l5], is the product of the number of phases and clusters’ diameter. 
Hence the approach of [39| appears to be inherently incapable to give rise to a constant time algorithm. 

Our deterministic variant of the network-decomposition procedure is the basis for our deterministic 
approximation algorithm for MDS. Our deterministic variant is closer to the algorithm of |3] than our 
randomized one. The main difference between our deterministic variant and the algorithm of [3| is that 
we use a different much faster coloring procedure for the supergraph induced by the ultimate partition. 
1.6 Related Work 

Network-decompositions for general graphs were studied in mmm- Dubhashi et al. m used network 
decompositions for constructing low-stretch dominating sets. Recently, Kutten et al. [35| extended Linial- 
Saks network-decompositions to hypergraphs. Many authors [25l [33l SH] studied network-decompositions 
for graphs with bounded growth. Distributed approximation algorithms is a vivid research area. See, 
e.g., [42] and the references therein. Distributed graph coloring is also a very active research area. 
See a recent monograph |9|, and the references therein. Schneider et al. m devised a distributed 
coloring algorithm whose performance depends on the chromatic number of the input graph. However, 
the algorithm of [U] provides no non-trivial approximation guarantee. To the best of our knowledge 
there are no known distributed approximation algorithms for the minimum t-spanner problem. Efficient 
distributed algorithms for constructing sparse undirected spanners can be found in |2n[ll7| . For centralized 
approximation algorithms for the minimum t-spanner problem, see [3011221 dl]. 

2 Preliminaries 

For a subset V C V, the graph G{V') denotes the subgraph of G induced by V. The degree of a vertex 
u in a graph G = {V,E), denoted degg(u), is the number of edges incident on v. A vertex u such that 
{u,v) G E is called a neighbor of v in G. The neighborhood of v in G, denoted rG'(u), is the set of 
neighbors of v in G. If the graph G can be understood from context, then we omit the underscript g- 
For a vertex v G V, the set v U r(I/) is denoted by r+(u). For a set VT C V, we denote by r'''(IT) the 
set W U U«;eiv r(r(;). The distance between a pair of vertices u,v G V, denoted distG'(rt, v), is the length 
of the shortest path between u and v in G. The diameter of G is the maximum distance between a pair 
of vertices in G. The chromatic number x{G) of a graph G is the minimnm nnmber of colors that can be 
used in a proper coloring of the vertices of G. 

3 Network Decomposition 

3.1 Procedure Decompose 

In this section we devise an algorithm for computing an (0(1), 0(n’^))-network-decomposition in 0(1) 
rounds, for an arbitrarily small constant e > 0. More generally, our algorithm computes a (3^,0(/c • 
7 ^ 2 /fc . log^ rr))-network-decomposition Q in 0(3^ • log* n) rounds, for any positive parameter k,l < k < 
logn, along with an 0{k ■ • log^ n)-coloring (p of the supergraph induced by Q. (The log* n term 

can be eliminated from the running time at the expense of increasing the number of colors used by p 
by a multiplicative factor of log^*^ n, for an arbitrarily large constant t. We will later show that the 
multiplicative factor of k in the second parameter of the network decomposition can also be eliminated 
without affecting other parameters.) The algorithm is called Procedure Decompose. The procednre rnns 
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on some supergraph G = {V,E) of the original graph G. Each vertex C G E is a cluster (i.e., a subset 
of vertices) of the original graph G = {V,E), and different clusters are disjoint. Observe that generally 
it may happen that V ^ The procedure accepts as input the supergraph G, the number of 

vertices n of G, the parameter k, and an upper bound s on the number of vertices of the supergraph G. 
It also accepts as input two numerical parameters e and t. The parameter e > 0 is a sufficiently small 
positive constant and t > 0 is a sufficiently large integer constant. Initially the supergraph is G itself, 
with each vertex v forming a singleton cluster {u}. Hence initially it holds that n = s. The procedure 
is invoked recursively. After each invocation the current supergraph G is replaced with a supergraph on 
fewer vertices, and s is updated accordingly. The parameter n, however, remains unchanged throughout 
the entire execution.) As a result of an execution of Procedure Decompose every vertex u in G is assigned 
a label label{v). The value of label{v) is equal to the color ip{Gy) of the cluster G^ of Q which contains v. 

Procedure Decompose partitions the graph G into two vertex-disjoint subgraphs with certain helpful 
properties. Specifically, one of the subgraphs has a sufficiently small maximum degree that allows us to 
compute a network decomposition in it directly and efficiently. The other subgraph can be partitioned into 
a sufficiently small number of clusters with bounded diameter. The latter property is used to construct 
a supergraph whose vertices are formed from the clusters. Since the number of clusters is sufficiently 
small, the number of vertices of the supergraph is small as well. Then our algorithm proceeds recursively 
to compute a network decomposition of the new supergraph, using fresh labels that have not been used 
yet. The recursion continues for k levels. Then each vertex is assigned the label of the supernode it 
belongs to. (Supernodes of distinct recursion levels may be nested one inside the other. In this case an 
inner supernode receives the label of an outer supernode. A vertex of the original graph G receives the 
(same) label of all supernodes it belongs to. Notice that a vertex belongs to exactly one supernode in 
each recursion level.) This completes the description of the algorithm. Its pseudocode is provided below. 
(See Algorithm [U) 

The algorithm employs two auxiliary procedures that we describe in detail in Section 13.21 The 
procedures succeed with high probability, i.e., with probability 1 — l/u'^, for an arbitrarily large constant 
c. The first procedure is called Procedure Dec-Small. It accepts a graph G with at most n vertices and 
maximum degree at most d. Procedure Dec-Small accepts also as input two numerical parameters, e and 
t, which are relayed to it from Procedure Decompose. Recall that e > 0 is a sufficiently small constant and 
t is a sufficiently large integer constant. The procedure computes an 0(min{d • (i^})-coloring of G in 
0(log* n) time. (The time is 0(1) if d > n^. Another variant of this procedure computes an 0{d? log^*^ n)- 
coloring in 0{t) time, for an arbitrarily large positive integer t.) Observe that for any integer p > 0, a 
proper p-coloring of a graph G is also a (O,p)-network-decomposition of G. (There are p labels, and each 
cluster consists of a single vertex. Thus the diameter of the decomposition is 0.) Procedure Dec-Small 
returns a (O,p)-network-decomposition S on line 5. It also returns a labeling function labels for vertices 
of a subset A. (We will soon describe how this subset is obtained.) The labeling labels also serves as a 
proper coloring for the supergraph induced by S. 

The second procedure which is invoked by our algorithm is called Procedure Partition. This random¬ 
ized procedure accepts as input an s-vertex supergraph G = {V,E) and a parameter q < 2 l\ogn ' 
partitions V into two subsets A and B, such that G{A) and G{B) have the following properties. The 
subgraph G{A) has maximum degree 0{q\ogn). The subgraph G{B) consists of 0{\V\/q) = 0{s/q) 
clusters of diameter at most 2 with respect to G. The procedure contracts each such cluster into a su¬ 
pernode. Let B denote the resulting set of supernodes and G{13) = {B,£{B)) the resulting supergraph. 
Specifically, the vertex set of Q{B) is H, and its edge set is £{B) = {(G, G') | G, G' G B, Bu G G,u' & 
C', such that {u,u') G E}. Procedure Partition returns the subset ACV and the set of supernodes B. 

The clusters in B are obtained by computing a dominating set D of R of size 0{\V\/q). Each vertex 
in D becomes a leader of a distinct cluster. Each vertex m B\D selects an arbitrary neighbor in D and 
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joins the cluster of this neighbor. Consequently, in all clusters all vertices are at distance at most 1 from 
the leader of their cluster. Hence all clusters have diameter at most 2. Initially, each vertex of V joins the 
set D with probability 1/q. Then the set B is formed by the vertices of D and their neighbors. Finally, 
the set A is formed by the remaining vertices, i.e., A = V \ B. In this stage the procedure returns the 
set of nodes A and the set of supernodes B which is obtained from B, and terminates. This completes 
the description of Procedure Partition. 


Algorithm 1 Procedure Decompose(G, n, k, s, e, t) 

1 : if s < 2c ■ logn then 
2 : return Dec-Smal^G, n, s, e, t) 

/* Compute directly a (0, 0(s^))-network-decomposition of G. (See Section [321) */ 

3: else 

4: {A, B) := Partition(G', q := v}/^) 

/* Partition G into A and B. (See Section [3.21 ) The maximum degree of G{A) is 0(n^/^ log n).*/ 
5: {S, labels) ■= Dec-Small(G(A), n, log n, e, t) 

/* Compute directly a ( 0 , 0(n^/^ • log^ n))-network-decomposition of G{A). (See Section [321) */ 
6 : (L, labels) ■= Decompose(^(H), n, k, ^j^) 

/* A recursive invocation on the supergraph G{B) that contains at most supernodes. */ 

7: for each vertex v of G, in parallel, do 

8 : if u E S' then 

9: label{v) := labels{v) 

10 : else ii V £ L then 

11 : label{v) := labelsiv) + A 

/* A = 7 • • log^ nj, where 7 is a sufficiently large constant to be determined later. */ 

12 : end if 

/* The labeling function label on S U L is defined by: for a cluster G £ S (respectively, G £ L) 
it applies to it the function labelsO (resp., labelsO + A). */ 

13: end for 

14: return (S U L, label) 

15: end if 


The recursive invocation of Procedure Decompose on line 6 returns a network decomposition L for 
the supergraph G{B). The for-loop (lines 7-13) adds (in parallel) A = 7 • log^ nJ to the color of each 
cluster of the network decomposition Lq of Q{B), where 7 is a sufficiently large constant to be determined 
later. Since the number of colors used in each recursive level is at most A, this loop guarantees that colors 
used for clusters created on different recursion levels are different. This is because the labeling returned 
by procedure Dec-Small on line 5 for clusters of S employs the palette [A] while the labeling computed 
in lines 10 - 12 for clusters of L employs labels which are greater than A. The termination condition of 
the procedure is the case s = 0(n^/^logn), i.e., when the number s of vertices in the supergraph G is 
already small. At this point the maximum degree of G is small as well (at most s — 1), and so coloring 
the supergraph (by Procedure Dec-Small) results in a sufficiently good network decomposition. 

Observe that our main algorithm will invoke the procedure on the original graph G. Hence in the 
hrst level of the recursion G = G, and each supernode is actually a node of G. In the second recursion 
level it is executed on the supernodes of nodes of the original graph G. In the third level it is executed on 
supernodes of supernodes, etc. Consequently, starting from the second recursion level supernodes have 
to be simulated using original nodes of the network. To this end each cluster that forms a supernode 
selects a leader which is used for simulating the supernodes. Moreover, the leader is used to simulate all 
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nested supernodes to which it belongs. Our supernodes are obtained by at most k levels of nesting. In 
each level of nesting a supernode is a cluster of diameter at most 2 in a graph whose nodes are lower-level 
supernodes. Hence a simulation of a single round on such a supergraph will require up to 3^^^ rounds. 

Next we provide several lemmas that will be used for the analysis of the algorithm. We leave the 
parameters e and t unspecified in all lemmas in this section, because they have no effect on the analysis. 

Lemma 3.1. Consider an invocation of Procedure Decompose on the original graph G = {y,E) with 
parameters n = \ V\, k, and s = n, for 1 < k < logn. The number of recursion levels in the execution of 
this Procedure (i.e., Decompose(G := G,n,k,s ■=n)) is k. 

Proof. In recursion level i, i = 1, 2,..., the parameter s is equal to Hence, in recursion level k 

the parameter s is equal to , and the recursion reaches the termination condition. (See lines 1-2 of 
Algorithm [U) □ 

Lemma 3.2. The number of labels used in the invocation of Procedure Decompose(G := G,n,k,s := n) 
is 0{k ■ rxf!^ log^ n). 

Proof. We show that the number of labels is j ■ k ■ log^ nj, where 7 is a sufficiently large constant. 
Specifically, the constant 7 needs to be larger than the constants hidden by the 0-notation in comments 
on lines 2 and 5 of the algorithm. (Recall that line 2 computes a (0, 0(s^))-network-decomposition, 
and line 5 computes a (0, 0(n^/^ log^ n))-network-decomposition. The constant 7 appears in line 11 of 
Algorithm dJ) The proof is by induction on ii = k — i + 1, where i is the recursion level. In other words, 
this is an inverse recursion on the number of recursion levels. For each index i G [fc], denote by Gj the 
supergraph on which Procedure Decompose is invoked on the fth level of the recursion. Note that at this 
point s = The inductive claim is that the zth level invocation of Procedure Decompose (on 

the supergraph Gj) employs at most • 7 • log^ nJ < (A: — i -|- 1 ) • 7 • log^ n labels. 

Base {ii = 1, i.e., i = k): In this case s < 2c • n^^^logn, the termination condition of the recursion 
holds, and thus the number of labels used in the decomposition is O(s^). (See line 2 of Algorithm [TJ) By 
the choice of 7 , the number of labels is at most 7 • ■ log^ n. 

Step: Suppose that the invocation has returned from level i -|- 1 of the recursion, and it is now at level 
i. By the induction hypothesis, line 6 of Algorithm [T] returns a labeling with 7 • ii^i ■ log^ nJ = 
7 • (£j — 1 ) • log^ nJ labels. Once line 11 is executed, the number of labels becomes 7 Tj • log^ nJ. 
This proves the inductive claim. 

In the end of recursion level i = 1 the algorithm terminates (after returning from all recursive invo¬ 
cations). In this stage it holds that ii = k, and the claim follows. □ 

Lemma 3.3. Each cluster created by the invocation above has diameter at most 3^“^ — 1. 

Proof. We prove by induction on.^ = A: — where i is the recursion level, that level-i clusters have 

diameter at most 3^“^ — 1 . 

Base {i = 1, i.e., i = k): In this case a (0, 0(s^))-network-decomposition is computed directly, and thus 
the diameter of all clusters in the graph G on which it is executed is 0 . (Recall that the argument G in 
the level-A; invocation is a supergraph of the original input graph G.) 

Step: First, observe that a (0,O(n^/^ • log^ n))-network-decomposition of G{A) is computed directly in 
line 5 of Algorithm [TJ Hence S consists of clusters of diameter 0 (with respect to supernodes of the 
supergraph G of the current recursion level). Next, we analyze the diameter of clusters in L. By the 
induction hypothesis, line 6 of Algorithm [T] (i.e., the recursive invocation of Procedure Decompose) returns 
a network decomposition in which all clusters have diameter at most 3^“^ — 1. This is a decomposition 
of the supergraph G{B). Consider a cluster C of diameter at most 3^“^ — 1 in G{B). Let x,y he a pair 
of vertices of G that belong to supernodes in C. Let G(x),G(y) G C be these two supernodes (clusters). 


such that X G C'(x), y G C{y). Since the diameter of C in Q{B) is at most 3^“^ — 1, there exist clusters 
C{x) = Ci,C 2 , = C{y) G C, such that t < 3^“^ — 1, and the following holds. There exist edges 

Cl = {ui,vi),e 2 = {u 2 ,V 2 ), = {ut-i,vt-i), such that for every i £ [t — 2],Ui G Cj, and for 

every i G [2, t — 1], Uj G Q+i. (See Figured] for an illustration.) By construction, each of the clusters 
Cl, C 2 ,..., Ct has diameter at most 2. Hence for i G [t — 2], it holds that distg(uj, < 2. Therefore, 

distg(x,y) < dist(^(x, til) + 1 + distg(ui, U 2 ) + 1 + distg(u 2 ) 1 * 3 ) + 1 + 

... +dist(5,(t>t_2,Mt-i) + 1 + dist(5(ut_i,y) < 2-t + t-l = 3-t-l. 

Since t < 3^“^, it follows that dist(j(x,y) < 3 • 3^“^ — 1 = 3^“^ — 1. Therefore, the diameter of C in G is 
at most 3^“^ — 1. Hence all clusters in S' U L have diameter at most 3^“^ — 1 in G. Since the number of 
recursion levels is k, the claim follows. □ 



Fig. 1. The clusters C{x) = Cl, C 2 , = C{y). 

Lemma 3.4. Suppose that all invocations of auxiliary procedures of Procedure Decompose have succeeded. 
Then the invocation computes a (3^“^ — l,0{k ■ ■ \o^ n))-network-decomposition. 

Proof. Consider a pair of distinct adjacent clusters C,C' G S U L. If C G S and C £ L then label{C) = 
labelsiC) £ [A], while label{C') = labehiC) > A. Hence label{C) + label{C'). 

If C,C' £ S then since Procedure Dec-Small returns on line 5 a network decomposition with a proper 
labeling labels{-), it follows that labels{C) / labels{C'), and so label{C) / label{C'). 

Finally, if C,C' £ L then inductively we conclude that labelL{C) / labeliiC'), and thus label{C) / 
label{C') too. (The induction base is the recursion level k, where the correctness follows from the 
correctness of Procedure Dec-Small invoked on line 2 of Algorithm dl) 
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Hence Procedure Decompose returns a partition S' U L into clusters of diameter at most 3^“^ — 1 (by 
Lemma l3.3p . and a proper labeling label{-) of this partition. By Lemma 13.21 the number of labels used 
by the labeling label{-) is 0{k ■ ■ log^n). Hence S U L is a (3^“^ — l,0{k ■ ■ log^ n))-network- 

decomposition for G, and label{-) is a proper labeling for the network decomposition S U L. □ 

Recall that the auxiliary procedures Dec-Small and Partition succeed with probability 1 — 1/n^, for 
an arbitrarily large constant c. Each of these procedures is invoked at most k < logn times during the 
execution of Procedure Decompose. Therefore, the probability that all executions of Procedure Dec-Small 
and Procedure Partition succeed is at least (1 — ~ 1 — Since c is an arbitrarily large 

constant, all executions of the auxiliary procedures succeed, with high probability. Hence Procedure 
Decompose computes a (3*^, 0{k ■ ■ log^ n))-network-decomposition, with high probability. 

The next lemma analyzes the running time of the algorithm. 

Lemma 3.5. Let Tpart{'n,q) (respectively, Tdecin,d)) denote the running time of Procedure Partition 
invoked with parameters n and q (resp., Procedure Dec-Small invoked with parameters n and d). We will 
assume that both these running times are monotone non-decreasing in both parameters. Then the running 
time of Procedure Decompose is 0{3^ ■ {Tpartin, -|- Tuedn, 2c ■ log n))). 

Proof. During the execution of Procedure Decompose the Procedure Dec-Small is executed k times, and 
Procedure Partition is executed k — 1 times. For i = 1,2, ...,k — l, in recursion level i both procedures are 
executed on supergraphs whose supernodes constitute subgraphs of diameter at most 3* of the original 
graph. Thus, the number of rounds required in level i is the product of the number of steps required to 
execute the procedure on the supergraph and the maximum diameter of supernodes. This running time 
is at most {Tpartin, + Tdec{n,n^^^^ogn) -\- 0(1)) • 3*. The running time of the last recursion level 
k in which the termination condition holds is Tdec{n,2c ■ n^/^logn) ■ 3^. Therefore, the overall running 
time is 

0(Eti3' • {Tpart{n,ny^) + Tdec{2c-n^lHogn))) = 0(3^ • {Tpart{n,n^/'^) + Tdec{n,2c ■ n^/'^Xogn))). □ 

Procedure Dec-Small and Procedure Partition are provided and analyzed in Section 13.21 Next we 
state the main results obtained by plugging these procedures into Procedure Decompose. See Section [32] 
for the proofs. 

Theorem 3.6. For any parameter k,l < k < logn, Procedure Decompose computes a {3^,0{k ■ ■ 

log^ n))-network-decomposition along with the corresponding 0{k ■ rifl^ ■ log^ n)-labeling function in time 
0(3^ • log*n), with high probability. Alternatively, one can also have the second parameter equal to 
0{k ■ r?!^ logn) and the running time 0{3^ ■ k). 

It follows that, an (0(1), n'^)-network-decomposition of an arbitrary n-vertex graph along with a 
proper n'^-labeling for it can be computed by a randomized algorithm, in 0(1) time, with high probability. 
See Section m 

3.2 Procedure Dec-Small and Procedure Partition 

We start with the description of Procedure Dec-Small. This procedure accepts a graph G with at most 
n vertices and maximum degree at most d, and computes an 0(min{d • n^, (i^})-coloring of G, where 
e is a fixed arbitrarily small positive constant. In other words, if d < n^ then an 0((i^)-coloring is 
computed, and otherwise an 0{d- n^)-coloring is computed. For computing an 0((i^)-coloring, Procedure 
Dec-Small employs the deterministic algorithm of Linial |38] that computes an 0(A^)-coloring of graphs 
with maximum degree A within 0(log* n) time. For computing an 0{d ■ n'^)-coloring. Procedure Dec- 
Small employs the randomized algorithm of Barenboim [5] that computes, with high probability, an 
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0(A • n'^)-coloring in 0(l/e) time, for an arbitrarily small e > 0. We henceforth refer to this algorithm 
as Procedure Random-Color. This completes the description of Procedure Dec-Small. Its pseudocode is 
provided below. 

For completeness, we provide a high-level description of the algorithm of Linial [38] and the algorithm 
of Barenboim |5] . The algorithm of Linial |38j starts with a legal n-coloring of the input graph obtained 
from the IDs of the vertices. It proceeds in phases, each of which reduces the number of colors while 
preserving the legality of the coloring. In each round the number of colors is reduced from p to 0(A^ logp), 
where p is the number of colors in the beginning of a round. (Initially p = n.) In the last round the 
number of colors is reduced from 0(A • polylog(A)) to O(A^). Each phase requires just a single round, 
and the overall running time of the algorithm of Linial is 0(log* n). (It is actually log* n + 0(1), but this 
precision is immaterial for our purposes.) 

Observe also that one can run Linial’s algorithm for just t rounds, for some positive integer parameter 
t, and obtain an 0(A^ Tog^*^ n)-coloring. For a single phase of Linial’s algorithm it employs A-union free 
set systems from the paper by Erdos, Frankel and Furedi [23]. A family F of sets over a given ground-set 
X is said to be A-union-free if for every A -|- 1 sets So, ^i,..., Sa G F, it holds that Sq ^ Erdos 

et al. showed that for any positive integers p and A, p > A -|- 1, there exists a A-union-free family F of 
p subsets over a ground-set X of size |Ai| = 0(A^ logn). 

Let be a proper p-coloring of G in the beginning of a phase of Linial’s algorithm. The algorithm 
associates a set Sc from F with each color c of p. Every vertex v that runs the algorithm computes a 
new color d G \ Such a color exists since T" is a A-union-free family, and (p(u) 7 ^ p{v) 

for every u G F(u).) The vertex v sets its new color <p'(u) by <p'(u) = d . Since d G X, \X\ = O(A^logn), 
it follows that is an 0(A^ logn)-coloring. Also, consider a pair of neighbors v and u. Observe that 
ip'iv) G \ while ^p'{u) G Hence ip'{v) 7 ^ p'{u), and thus ip' is a proper coloring. See [38] 

or [9] Chapter 3.10 for more details. 

The algorithm of Barenboim [5] (Procedure Random-Color) proceeds in phases as well, however, each 
phase consists of a randomized procedure. In this procedure each vertex either succeeds in selecting a 
final color, or fails and continues to the next phase. Each vertex selects a color from the range [A • n^] 
uniformly at random, where e > 0 is an arbitrarily small constant. A vertex succeeds if and only if it 
selects a color that has not been selected by any of its neighbors (either in the current round, or as a final 
color in a previous round). Otherwise, it fails, discards its color, and continues to the next phase. Hence, 
the probability that a vertex v fails to select a color that is different from the colors of its neighbors 
is at most l/n*^. If we run this procedure for \c/e\ rounds, for a sufficiently large constant c, then the 
probability that a given vertex v fails on all these rounds is at most 1/n'’. Hence, by union bound, after 
[c/e] rounds all vertices succeed with probability at least 1 — i.e., with high probability. 


Algorithm 2 Procedure Dec-Small(G, n, d, e, t) 

1 : if d < then 

2 : compute an 0((i^)-coloring of G using the algorithm of Linial [38] 

/* alternatively, one can compute here an 0(d^log^*^ n)-coloring in 0{t) time */ 
3: else 

4: compute an 0{d ■ n'^)-coloring of G using Procedure Random-Color 

5: end if 


Observe that if < d ■ then d < n’' and an 0(d^)-coloring is computed. Otherwise, dd > d ■ n^, 
and an 0{d ■ n^)-coloring is computed. Therefore, the number of colors is 0(min{d^,d • n^}). Recall also 
that the running time of the algorithm of Linial [38] is 0{log* n), and the running time of Procedure 
Random-Color is 0(1). Thus we obtain the following lemma. 
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Lemma 3.7. Procedure Dec-Small invoked on a graph G with maximum degree at most d computes, with 
high probability, an 0 (min{d^, d ■ coloring, which is a ( 0 , 0{m.m.{d^, d ■ n'^}))-network-decomposition. 
If d < n^, the running time of Procedure Dec-Small is 0(log*re). Otherwise, it is 0(1). 

Another variant of the Procedure Dec-Small checks if d < ^ , and if it is the case it invokes the 

t-round version of Linial’s algorithm. Otherwise it invokes line 4 of Algorithmic] (i.e., the algorithm from 
[5]). This modified procedure always requires constant time. (Assuming that t = 0(1).) If d < ^ ; 

it computes an 0{d? log^*^ n)-coloring. Otherwise (in this case d^ log^*^ n > d ■ n'^) it computes a (d • n^)- 
coloring. To summarize: 

Lemma 3.8. A modified variant of Procedure Dec-Small computes an 0(min{d^ log^*^ n,d-n‘^})-coloring 
in 0{t) time, with high probability. In particular, the running time is constant if t = 0(1). 

Next, we describe Procedure Partition. Procedure Partition accepts as input a graph G = {V, E) and 
a positive parameter q, and partitions V into two subsets A and B, such that G{A) and G{B) have the 
following properties. The subgraph G{A) has maximum degree O(glogn). The subgraph G{B) consists 
of 0{\V\/q) clusters of diameter at most 2. The procedure contracts the clusters of B into supernodes, 
which form the supergraph Q{B). The clusters in B are obtained by computing a dominating set D of B 
of size 0{\V\/q). Each vertex in D becomes a leader of a distinct cluster. Each vertex in B\D selects an 
arbitrary neighbor in D and joins the cluster of this neighbor. Consequently, in all clusters all vertices 
are at distance at most 1 from the leader of their cluster. Hence all clusters have diameter at most 2. 

Initially, each vertex of V joins the set D with probability 1/q. Then the set B is formed by the 
vertices of D and their neighbors. Finally, the set A is formed by the remaining vertices, i.e., A = V \ B. 
In this stage the procedure returns the set of nodes A and the set of supernodes B which is obtained from 
B, and terminates. This completes the description of the procedure. Its pseudocode is provided below. 


Algorithm 3 Procedure Partition(G, q) 

An algorithm for each vertex u G P. 

1: V joins D with probability 1/q and informs its neighbors 
2 : if V has joined D or a neighbor of v has joined D then 
3: V joins B 

4: else 

5: V joins A 

6 : end if 

7: if u € D then 

8: V initializes a singleton cluster C^, becomes the leader of G^, and sends Id{v) to all neighbors in B 

9: end if 

10: if u G -B and v receives at least one message from a leader of a cluster in D then 
11: V joins a cluster of an arbitrary neighbor in D 

12 : end if 

13: B ■.= the set of supernodes obtained by contracting all clusters G^ C B 
14: return (A,B) 


Note that the sets A and B are returned in a distributed manner. In other words, each vertex knows 
whether it belongs to A or to B. If it belongs to B, then it knows the identity of the leader of its cluster. 
The leaders of the clusters represent the supernodes formed by the clusters. Thus, when we say that 
a supernode performs some action, it is actually performed by the leader of the cluster that forms the 
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supernode. (For nested supernodes the operations are performed by the leaders of the innermost clusters, 
which are vertices in the original input graph G.) 

In the next lemmas we prove that Algorithm [3] computes a partition with the properties described 
above. 

Lemma 3.9. Suppose that Procedure Partition is invoked on a graph G = (F, E) and a positive parameter 
q. Then the subgraph G{A) induced by the set A which is returned by the procedure has maximum degree 
0{q ■ logn), with high probability. 

Proof. Consider a vertex v £V such that v has at least c - q - Inn neighbors in G, for a sufficiently large 
constant c. Denote 5 = degQ(n). Let yi,y 2 , .■.,y 5 denote the neighbors of v in G, and let yo = v. The 
probability that none of these neighbors join D is 

Pr(yi i D, ...,ys i D) = Pr(y, ^ D) = nti(l - l/q) = (1 - l/qf < (1 - < l/nC 

Hence, by union bound, the probability that at least one vertex v with at least c - q - Inn neighbors does 
not have a neighbor in D is at most Hence with probability at least 1 — all high-degree 

vertices (vertices with degQ(n) > c ■ q ■ Inn) end up in B. Hence, with high probability, the maximum 
degree of a vertex in A = F \ H is 0{q\ogn). □ 

The next lemma analyzes the number of supernodes and their diameters. 

Lemma 3.10. Suppose that Procedure Partition is invoked on a graph G = {V, E) and a parameter 
q < 2 cAogn ’ some constant c > 1. Then the set B returned by the procedure has the following 
properties. With high probability, B consists of 0{\V\/q) supernodes. All supernodes of B are clusters of 
diameter at most 2 in G. 

Proof. Recall that the set B is created by contracting the clusters of B into supernodes. First, we prove 
that all clusters of B have diameter at most 2. Let C be a cluster of B. Let u,v G C be any pair of 
vertices in the cluster. Then either one of these vertices is the leader of the cluster and dist(ri, u) = 1 
or both u and v are connected to the same leader, and so dist(M, v) < 2. (Since u and v belong to B 
they must have a leader neighbor, unless they are leaders themselves. Since u and v belong to the same 
cluster, and there is exactly one leader in each cluster, if u and v are not the leaders, they are connected 
to the same leader.) Next, we prove that B consists of 0{\V\/q) supernodes. Note that the number of 
supernodes in B is equal to the number of vertices in D, since each vertex in D becomes a leader of a 
cluster (i.e., of a supernode). Let X denote a random variable that counts the number of vertices in 
D. Since each vertex in F joins D with probability l/q independently of other vertices, it holds that 
1E(A) = |F|/g. By Chernoff bound for upper tails (see, e.g., |H], Chapter 4), 

Pr[X > 2\V\/q] = Pr[X > 21E{X)] < (e/4)®^(^) = (e/4)l^l/'? < (e/4)2'=-'°g^ < l/nC □ 

Finally, note that each line of Procedure Partition is either performed locally, or involves sending 
messages to neighbors. The latter requires one time unit. Therefore, the running time of Procedure 
Partition is 0(1). 

Lemma 3.11. Procedure Partition requires 0(1) time. 

Combining Lemmas 13.31 - [331 with Lemmas 13.71 - [3T0] imply the following results. 

Theorem 3.12. For any parameter k,l < k < logn, Procedure Decompose computes a {‘i^,0{k ■ ■ 

log^ n))-network-decomposition along with the corresponding 0{k ■ r?!^ ■ \og^ n)-labeling function in time 
0(3^ • log* n), with high probability. 
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Consider now a variant of Algorithm [T] (Procedure Decompose) in which in Procedure Dec-Small we 
always invoke Procedure Random-Color with a parameter e. (As opposed to Algorithm [2] where we do it 
only when d < n^.) Also, in Algorithm [ 1 ] we now set A •(— (2 • c • ■ logn) • n^. 

Then, by the previous argument, this modified variant of Procedure Decompose computes a {3^,0{k ■ 
^i/k+e jQg 77 ,)).network-decomposition along with a legal log n)-labeling function in time 0(3^/e). 

By substituting e = 1/k we conclude: 

Theorem 3.13. A ,0{k ■ n^^^logn))-network-decomposition along with the appropriate proper (with 
respect to this decomposition) 0{k ■ logn)-labeling can he computed in 0{k ■ 3^) time, with high 

probability. 

In particular, by setting k to be an arbitrarily large constant we obtain an (0(1), n'^)-network- 
decomposition along with a proper n'^-coloring for it in randomized constant time, for an arbitrarily 
small constant <5 > 0. 

Corollary 3.14. An {0{l),n^)-network-decomposition of an arbitrary n-vertex graph along with a proper 
-labeling for it can he computed by a randomized algorithm, in 0 ( 1 ) time, with high probability. 

A yet another variant of Procedure Decompose (Algorithm [T]) is obtained if in Proc Dec-Small we 
always invoke the t-round variant of Linial’s algorithm [38], for some positive integer parameter t. (Again 
we do it regardless of the value of d.) Also, for this variant we set A = 7 • log^ n log^*^ re, where 
7 > 0 is a sufficiently large constant. (Specifically, the t-round variant of Linial’s algorithm computes an 
0(A^ log^*^ re)-coloring of the input re-vertex graph with maximum degree A. The constant 7 should be 
larger than the constant hidden by the O-notation in 0 (Alog^*^ re).) 

The resulting algorithm computes a (3^,0(fc • re^/^ • log^ re • log^*^ re))-network-decomposition with a 
proper 0{k ■ nf!^ ■ log^ re • log^*^)-labehng for it, in 0(3^ • t) time. 

Corollary 3.15. For any n-vertex graph G and parameters k = 1,2,...; t = 1,2,...; e > 0, one can 
compute a (3*’, 0{k ■ ■ log n))-network-decomposition (respectively, (3^, 0{k ■ nf!^ ■ log^ re • log^*^ re))- 

network-decomposition) with an appropriate labeling function in 0{3^/e) (resp., 0(3^ • t)) randomized 
time, 

4 Refining the Algorithm 

In this section we argue that one can save a factor of k in the number of labels, and compute a 
(3^, 0 (re^/^"*“'^ log re))-network-decomposition (and a (3^, 0(r?l^ log re log^*^ re))-network-decomposition) with 
an appropriate labeling in 0{‘3f je) (resp., in 0(3^ • t)) time. While this improvement is negligible when 
k is small, it becomes significant when k is superconstant. We remark, however, that in the context of 
the current paper we are mainly interested in the regime of small k. 

To describe this improvement we need the notions of arboricity and //-partition. (We refer the reader 
to | 6 | and [9] for a more elaborate discussion on this topic.) 

The arboricity a{G) of a graph G = {V, E) is the minimum number t of edge-disjoint forests Fi, T 2 ,..., Ft, 
such that E = An H-partition of G = {V,E) with degree at most A, for some 

number A, is a partition of the vertex set V of G into vertex disjoint subsets V = Hi Ci Hj = $ 

for every pair of distinct indices i ^ j, i,j G [/], such that for every index i G [/] and every vertex v & Hi, 
the number of neighbors deg(u, U^^jFj) that v has in H-sets Hj with an index j > i is at most A. 

Consider again Procedure Decompose. (See Algorithm (U) For i = l,2,...,k, let Gi denote the 
supergraph on which the procedure is invoked in the ith. level of recursion. In particular, Gi = G is the 
original graph. Also, in all levels i < k — 1 the procedure enters lines 4 - 15, and in the last level i = k 


14 


it enters the termination condition (line 2). In the former case Procedure Decompose invokes Procedure 
Dec-Small (in line 5), which returns the collection S of clusters. (It also returns the labeling function that 
is immaterial for the current discussion.) For i = 1,2, ...,k — 1, let Si denote the set of clusters returned 
by Procedure Dec-Small on line 5 of the ith level recursive invocation of Procedure Decompose. Finally, 
in level k of the recursion Procedure Dec-Small is invoked in line 2. Denote by the decomposition that 
it returns. 

Lemma 4.1. is the network decomposition that Procedure Decompose returns (in line 14 of the 

first level recursive invocation). Moreover, for any index i, k>i>l, is the network decomposition 

that the ith level recursive invocation of Procedure Decompose returns. 

Proof. The proof is by induction on i. 

Base (i = k): In this case Procedure Decompose returns the output Sk of an invocation of Procedure 
Dec-Small (on line 2 of Algorithm [1]) . 

Step: Consider some 1 < i < k. The ith level recursive invocation returns S' U L in line 14. Recall that 
S = Si is a network decomposition for G{A) computed in line 5 of Algorithm [TJ (In all levels except the 
first one S is actually equals to A. In the first level S = {{u} | u € A}.) Also, L (computed by the (i-l-l)st 
level recursive invocation of Procedure Decompose; see line 6 of Algorithm [1]) is a network decomposition 
for G{B). By induction hypothesis the latter is Uj^i_^_iSj. Hence S U L = Uj^iSj, as required. □ 

In the next lemma we show that (Si, S 2 ,..., Sk) is an Fl-partition with relatively small degree of the 
supergraph ^(U^T^Si) induced by the network decomposition 

Lemma 4.2. (Si, S 2 ,..., Sfc) is an H-partition with degree 0(n^/^logn) of the supergraph 0(ujL;^Sj). 

Proof. Again, by an induction on z, /c > f > 1, we show that (S,, Sj+i,..., Sfc) is an i4-partition of 
Gii->i=iSi) with maximum degree 0 (n^/^logn). 

Base (z = k): In this case we need to show that the maximum degree in G{Sk) is 0(n^/^logn). By 
the termination condition of Algorithm [ 1 ] (line 1 ), |Sfc| = 0(n^/^ logn), and thus the same upper bound 
applies to its maximum degree. 

Step: For some z, 1 < z < fe, we argue that for any cluster C G Sj its degree in G(}Jj=iSj) is log n). 

By Lemma Fd.ll U^^jSj is the network decomposition for Gfilj^iSj) that the zth level recursive invocation of 
Procedure Decompose returns. By construction, Sj is the set of clusters with degree at most Oirfil^ log re) 
in the supergraph Gi. Since clusters of are obtained by merging clusters of Gi, it follows that 

the degree of G in \J^-^^Sj is no greater that its degree in Gi, i.e., at most 0(re^/^logre). □ 

To recap. Lemma 14.21 shows that in addition to computing a network decomposition Q of its input 
graph G, Procedure Decompose also computes a low-degree Fl-partition of the induced supergraph G{Q). 
(Here Q = ujL^^S'i, and the Ll-partition is {Si, S 2 ,..., Sk). The degree of the partition is 0(re^/^ log re).) 

For the variant of Procedure Decompose that we describe in this section we do not actually need to 
explicitly compute the labeling function during the execution of the procedure. As a result Procedure 
Dec-Small can be greatly simplified. Specifically, if it is invoked on a supergraph G = {V,E), then 
it returns V as its output partition. If it is invoked on a subgraph G{U) = {U,E{U)) of the original 
graph G, then it returns a partition of U into singleton clusters, i.e., {{re} | re G U}. Observe that in 
this simplified form Procedure Dec-Small requires 0(1) time. As a result the overall running time of 
Procedure Decompose becomes 0(3^) rather than 0(3^/e) or 0{3^ ■ t). 

Next, we utilize the Fl-partition {Si, S 2 ,..., Sk) of G{Q) for computing an 0(re^/^ log^ re log^*^ re)- 
coloring of G{Q) in 0(3^ • t) time, or alternatively, an 0(71^/^^*^ log re)-coloring of G{Q) in 0(3^/e) time. 
Such colorings can be viewed as labelings of the network decomposition Ui^^Si = Q. (For every vertex 
z; G P, its label will be equal to the color of the cluster Gy G Q that contains it.) 
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To simplify presentation, consider an n-vertex graph G' = {V', E') and an //-partition (5i, 82 -, ■■■, Sk) 
for G' with degree A = log re). We will argue that G' can be efficiently colored. To implement 

this coloring in a supergraph G{Q), we will need to multiply the running time by the maximum diameter 
of a cluster in Q, i.e., by 0(3^). We start with arguing that G' can be colored in 0(^^log^*^ re) colors 
in 0{t) time, by a deterministic algorithm. (This algorithm is closely related to Algorithm Arb-Linial 
from [6], based on Linial’s algorithm [38]. The current algorithm is however more general than Algorithm 
Arb-Linial.) The algorithm starts by orienting all edges {u,v) in the following way: let iu (respectively, 
iy) be the index of the set Si^ (resp., Si^) which contains u (resp., v). If iy < iy then the edge is oriented 
towards v. If the opposite holds than it is oriented towards u. If iy = iy then the edge is oriented towards 
the endpoint with a greater Id. 

Observe that under this orientation each vertex v has at most A outgoing edges incident on it. The 
opposite endpoints of these edges will be referred to as the parents of v. Let n(u) denote the set of parents 
of V. 

Let be a proper p-coloring of G'. We argue that a legal 0(A^ logp)-coloring tp' of G' can be 
computed within one single round. To this end we again employ an A-union-free family E of p sets (due 
to [23|, see also Section [32] of this paper). Each color class c of (/J is associated with a set Xy G E. A 
vertex v computes a color ip'{v) which belongs to \ Such a color necessarily exists, 

because E is an A-union-free family. Also, for an edge {v,u) G E', suppose without loss of generality 
that u G n(u). Then ip'{v) G X^^(^y■^ \ while p'{u) G and so p'{v) 7 ^ p'{u). By [23], a family 

E over a ground-set of size O(A^logp) exists (and can be efficiently constructed). Thus, p' is a proper 
0{A? logp)-coloring. By repeating this recoloring step for t times, we obtain an 0{A^ log^*^ re)-coloring in 
0{t) rounds. (We start with an initial re-coloring of G'. Specifically, each vertex uses its Id as its initial 
color.) 

Corollary 4.3. Are 0(re^/*^ log^ re log^*^ re)-co/onre 5 f of Q{Q) can he computed in 0(3^ • t) time, for any 
f = l,2,.... 

Observe that this argument shows in fact that the arboricity of G{Q) is 0(A) = OirA/^logn), and 
thus Q is a (3^“^ — 1, log re))-network-decomposition. We summarize these results in the following 

corollary. 

Corollary 4.4. Procedure Decompose, invoked on an n-vertex graph G = {V, E) with a parameter k = 
1 , 2 ,..., computes a (3^“^ — I,0{n^/^\ogn))-network-decomposition Q and an H-partition {Si,...,Sk) of 
degree 0(re^/^logre) and length k for G{Q) in 0(3^) randomized time, with high probability. Moreover, 
for a parameter t = 1,2,..., one can compute in 0(3^ • t) time an Oinf!^ log^ re log^*^ n)-labeling for Q. In 
particular, by setting t = log* re one can get here time 0(3^ log* re) and labeling with log^ re) labels. 

Note that the 0(A^log^*^ re)-coloring algorithm for G' that was described above does not require the 
fact that the //-partition { 81 , 82 , ■.., 8 k) has small number of sets. Next we show that this //-partition 
can be used in a more explicit way to compute an 0{A ■ re^)-coloring of G' in 0{k/e) time. 

First, every vertex v of 8 k tosses a color (p{v) uniformly at random from the palette [A • re^j. It checks 
if its color is different from the colors of all its neighbors in 8 k. If it is the case, it finalizes its color. 
Otherwise, it tosses its color from the same palette again. The process is repeated for [c/e] rounds, for 
a sufficiently large constant c. As we have already seen (see Lemma [3771 and the discussion preceding it), 
in 0(l/e) rounds we will obtain a legal 0(A • re^)-coloring (pk for 8 k, with high probability. (Recall that 
the maximum degree in 8 k is at most A.) Define also (pk = Pk- 

Suppose that we have already computed an 0(A-re'^)-coloring ipi for G'{L)j^^ 8 j), for some i, 2 < i < k. 
Next we show how to extend this coloring into an 0(A • re'^)-coloring (pi-i for G'{Uj^^_i 8 j). To this end 
every vertex v G 8 i-i tosses a color from [A • n^] uniformly at random, and checks if its color is different 
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from the colors (either tossed on this round, or finalized colors) of its neighbors in If it is 

different from them, then v finalizes its color. Otherwise, it continues to the next round. The entire 
process continues for 0 (l/e) rounds. 

The key observation required for the analysis is that v G Si-i has at most A neighbors in 
Thus, a legal {A ■ n'^)-coloring (pi-i for 5j_i will be computed, with high probability, within additional 
0(l/e) rounds. It is then combined in a trivial way with the coloring ipi for to obtain the {A ■ n*^)- 

coloring for 

Theorem 4.5. Given an H-partition {Si, S 2 , ■■■, S^) with degree at most A for an n-vertex graph G' = 
iy' ,E'), and a parameter e> 0, an {A - n^)-coloring of G' ean be computed in 0{k/e) rounds. 

By invoking this algorithm on the network decomposition Q we obtain; 

Corollary 4.6. Using a — \,0{yf^ \ogn))-network-decomposition Q of the input graph G and an 

H-partition {Si, S 2 , ■■■, Sk) for Q with degree 0(n^/^ log n), one can compute an log n)-labeling 

for Q within 0{3^ ■ k/e) randomized time. 

By substituting e = 1/k we get: 

Corollary 4.7. A {3^~^ — 1, 0{y/^ log n))-network-deeomposition Q along with an 0{y^^ log n)-labeling 
for it ean be computed in 0{3^ ■ randomized time. 

5 Decompositions with a smaller number of labels 

When k is small the logarithmic factor in the number of labels {Olpf!^ logn)) of the network decomposi¬ 
tion Q from Corollary 14.71 is almost negligible. However, for large k (e.g., k = H(logn)) this logarithmic 
factor becomes dominant. In this section we describe a modification of our algorithm that produces 
{exp{0{k)}, ))-network-decomposition in exp{0{k)"\ ■ log^^^ n time. (For graphs of girth at least 6 

the running time of this algorithm is even better, specifically exp{0{k)} ■ exp{O{ Vloglog n)}.) Observe 
that for k = H(loglog n), the overhead factor of log^'^^ n can be swallowed by the O-notation in exp{0{k)}. 
This version of our algorithm is closely related to the deterministic algorithm of Awerbuch et al. [3]; 
in fact, our algorithms in this section can be viewed as a randomized version of their algorithm. Their 
deterministic algorithm requires time (logn)®(^\ and so we essentially show here that their algorithm 
can be made faster by means of randomization. 

The difference between the new variant of our algorithm (which we introduce here; we will refer to 
it as Procedure RS-Decompose) and the original version of our algorithm (described in Section 13.ip is a 
different algorithm for Procedure Partition. (See Algorithm [3l) The new variant of Procedure Partition 
which we will next describe will be called Procedure RS-Partition. (RS stands for the acronym of ’’ruling 
set”.) 

In a graph G = iV, E) for a vertex set [/ C R and positive integer parameters r, <5 a subset IT C [7 is 
called an {r, 6)-ruling set for U if the following two conditions hold: 

(a) Every pair of distinct vertices w,w' G IT satisfy distG(ic, tc') > 5. 

(b) For every vertex u G U there exists a ’’ruling vertex” (also called ’’ruler”) tc G IT such that 
distG(tc, u) < r. 

Observe that an MIS is a (2, l)-ruling set. In the description of Procedure RS-Partition we will 
assume that we have an efficient distributed subroutine for computing (r, 5)-ruhng sets for r = 3 and 

6 = 0(1). We will later elaborate on this subroutine. Procedure RS-Partition starts with computing a 
(3, (I)-ruhng set W for the set 1/ = {u G T | deg(ri) > q} of high-degree vertices of G. (Recall that q is an 
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input parameter of Procedure RS-Partition.) Then every vertex w sends an exploration message to 
distance 5. Every vertex v that receives an exploration message from two distinct rulers w' ,w" E W 
assigns himself to the ruler w which is closer to it. (Ties are broken in an arbitrary but consistent manner 
by comparing rulers’ identities.) 

As a result of these explorations clusters {Cyj | rc E TT} are formed. Observe that these clusters all 
have strong radius at most <5, and that every u & U (i.e., every high-degree vertex) is assigned to some 
cluster. (This collection of clusters is often called a ruling forest. See, e.g., [3].) Procedure RS-Partition 
now forms the set B of supernodes by contracting these clusters exactly as in line 13 of Algorithm [3l 
Further, it creates the set A by setting A •(— R \ {Uuj&wCw)-, i-e., every vertex v which is not clustered is 
assigned to A. Observe that for every u E A, it holds that deg(u) < q. Finally, Procedure RS-Partition 
returns the pair {A,B), exactly as in line 14 of Algorithm [3l 

Lemma 5.1. Suppose that Procedure RS-Partition is invoked on a graph G = {V,E) and a positive 
parameter q. Suppose further that it uses a subroutine for computing a {3,5)-ruling set, for a positive 
integer parameter 6. The the subgraph G{A) has maximum degree smaller q. Moreover, B consists of at 
most \ V\/q supernodes, each of which is a cluster of strong diameter at most 25. 

Proof. All the assertions of the lemma were already argued in the preceding discussion, except for the 
claim that \B\ < |R|/g. We next show this claim. Recall that every supernode of B originated from a 
cluster Cu), w E W, where VF is a (3, (i)-ruling set for the set U of vertices with degree at least q. Hence 
for two distinct clusters Cu,,Gw' from the collection B = {G^ | u) E W}, it holds that distG(rc, u;') > 3, 
and deg(u)), deg(t(;') > q. All (immediate) neighbors of w (respectively, w') are assigned to the cluster 
Gw (resp., Guji)-, and these sets of neighbors are disjoint. Hence IC^I > q for every w E W, and 
\B\ < \V\/q. □ 

We now use Procedure RS-Partition instead of Procedure Partition within Procedure RS-Decompose. 
The diameter of clusters in the modified procedure becomes {25 + 1)^ instead of 3^, but the factor 
logn is shaved from the bound on arboricity. (This is because the bound on deg(A) for A returned by 
Procedure RS-Partition is q instead of 0{q ■ logn). Hence as a result we obtain a ((25 -|- 1)^, 0(n^/^))- 
network-decomposition Q and an H-partition (5i, S 2 , ..., Sk) of degree 0{n^^^) of length k for Q{Q). (See 
Corollary 14.41 for a comparison.) 

To analyze the running time we need to specify the black-box procedure for computing a (3,5)-ruling 
set W for the set U of high degree vertices. Barenboim et al. [To] (based on [3T] and [12]) showed that 
(2, 2)-ruling sets can be computed in 0(log^/^ A)+exp{0{ Vlog log n)} randomized time in general graphs, 
and that (2, 3)-ruling sets can be computed in graphs with girth at least 6 in just exp{0 (Vlog log n)} time. 
By running their routine in G^ we guarantee that any two distinct vertices w,w' E W are at distance at 
least 2 in G^, i.e., at distance at least 3 in G. On the other hand, the domination parameter grows by a 
factor of 2, i.e., we obtain a (3,4)-ruling set in 0(log^/^ A) -|- exp{0{ Vlog log n)} time in general graphs, 
and a (3,6)-ruling set in exp{0{ Vlog log re)} time in graphs of girth at least 6. Hence the running time of 
Procedure RS-Partition becomes now O(log^'^^n) for general graphs and exp{0{ \/log logn)} for graphs 
of girth at least 6, instead of the running time of 0(1) for Procedure Partition. (See Lemma [3.11l l Hence 
by Lemma 13.51 the overall running time of Procedure RS-Decompose becomes 0((25 -|- 1)^ • log^^^ n) for 
general graphs, and 0((25 -|- 1)^ • exp{0{y^log log n)} for graphs of girth at least 6. In the former case 
5 = 4, while in the latter it is 6. To summarize, we have proved the following theorem. 

Theorem 5.2. Procedure RS-Decompose invoked on an n-vertex graph G = {V,E) with a parameter 
k = 1 , 2 ,..., computes an {{0{1))^,n^/^)-network-decomposition Q and an H-partition {Si, S 2 ,..., Sk) of 
degree A = 0(n^/^) of length k for G{Q) in (0(1))^ log^^^ n randomized time for general graphs, and in 
(0(1))*’ • exp{0{ Vlog log n)} time in graphs of girth at least 6. 
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See Corollary 14.41 for the comparison between the result here and the result that we have for the 
original variant of our algorithm. 

Also in a way analogous to Corollary 14.71 Theorem 15.21 implies that we can also compute a labeling for 
the network-decomposition Q. The time required to compute an 0{A ■ n^/^)-coloring for G{Q) given an 
//-partition as above is, by Theorem 14.51 at most 0{Diam{Q) ■ /c^) = (0(1))^ = exp{0{k)}. The number 
of labels (colors) is 0{A ■ = 0( We summarize the properties of the network-decomposition 

Q in the next corollary. 

Corollary 5.3. An {exp{0{k)},n^^^)-network-decomposition Q along with an -labeling for it can 

be computed in exp{0{k)} n (respectively, exp{0(A:)} • O(vTogTog^)} ) randomized time in general 

graphs (resp., in graphs of girth at least 6). 

Observe that randomization was used by the modified variant of Procedure Decompose only for com¬ 
puting a ruling set and for computing the labeling. There is a deterministic algorithm for computing 
(3, 0(log n))-ruhng sets in O(logn) time due to [3]. If we plug it in the above algorithm the diame¬ 
ter of Q grows from (0(1))^ to 0(logn)^“^), and consequently, the running grows to 0((logn)^“^) as 
well. (The most time-consuming step involves computing a (3,0(log n))-ruling set in the last phase of 
the algorithm, i.e., in a supergraph in which each cluster has diameter (0(logThis requires 
(0(logn))^“^ time.) Hence we obtain the following result, which is a generalization of the network de¬ 
composition of [3]. (They arrived to the same result with k = -^/lognlog logn, i.e., they obtained an 
{exp{O (Vlog n log log n )}, exp{O (-y/log n log log n)})-network-decomposition.) 

Corollary 5.4. An (((0(log n)Y~^-network-decomposition Q along with an H-partition {Si, S 2 , ■■■, Sk) 
of degree A = 0{n^^^) of length k for Q{Q) can be computed in deterministic (0(log n))^~^ time in general 
graphs. 

This also gives rise to a construction of 0((log n)^“^)-spanner with edges, in deterministic 

0((logn)^“^) time, in the CONGEST model. This is achieved by adding one edge for every pair of 
adjacent clusters of the decomposition of Corollary 15.41 By setting k = ^ constant c > 1, 

one can get 0{n^^‘^) time and 0{n ■ log'^n) edges. In particular, this results in a sparse skeleton (with 
n-polylog(n) edges), in time 0{n^), for an arbitrarily small constant e > 0, in the deterministic CONGEST 
model. 

Using the //-decomposition of G{Q) from Corollarv l5.4l an O(A^) = 0(n^/^)-labehng for it (i.e., 0{A^)- 
coloring) for G{Q) can be computed by Algorithm Arb-Linial within additional 0{Diam{Q) ■ log* n) = 
(0(logn))^“^ • log* n deterministic time. (This is another point in which this deterministic routine is 
different from that of [3]. To compute the coloring Awerbuch et al. [3] used here 0{Diam{Q) - A-logn) = 
0{n^/^ \og^ n) time, but the number of colors was only 0{n^/^) instead of 0{n ^/^). Since we insist on 
having a deterministic polylogarithmic time, this modification is crucial.) 

Corollary 5.5. For any positive integer k, an {{0{\ogn))^~^,n^/^)-network-decomposition Q along with 
an 0{n^^^)-labeling of it can be computed in (0(logn))^“^ • log* n deterministic time. 

By running a /-round version of Algorithm Arb-Linial, for some positive integer constant t, one can 
also have here running time 0((logn)^“^ • t) = 0((log n)^“^), but the number of colors (labels) becomes 
0{r?l^ log(i) 

6 Applications 

6.1 An Approximation Algorithm for the Coloring Problem 

The results described in the previous sections (Theroem 13.61 see also Corollary I4.7jl imply an approx¬ 
imation algorithm for the optimization variant of the coloring problem. A distributed approximation 
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algorithm for the graph coloring problem (based on an ( 0 ( 1 ), 0 (n^/^'*'’^))-network decomposition) was 
given in [5]. We describe here a generalization of that algorithm which works with any network- 
decomposition. The algorithm starts by computing a (3^ — 1, 0(n^/^ log n))-network-decomposition Q 
with an log n)-labeling label{-) for it. See Corollarv l4.7l Then in each cluster C the entire induced 

subgraph G{C) is collected into the leader vertex vc of C. The leader vertex vc computes locally the 
optimum coloring ipc for C. Finally, vc broadcasts (a table representation of ipc) to all vertices of C. 

Each vertex u that receives this broadcast computes its final color ^^{u) by V’(^) = {'fc{u),label{u)). The 
running time of this algorithm is the sum of the time required to compute the decomposition Q (i.e., 

0(3^ • fe^)) with the time required for the computation of the colorings (pc- The latter is dominated by 
the diameter of Q, times a small constant. The overall running time is therefore 0(3^ • k"^). 

The next lemma shows that the coloring ^ provides an 0(n^/^ log n)-approximation to the optimal 
coloring for G. 

Lemma 6.1. ijj is a proper logn ■ x{G))-coloring. 

Proof. Consider an edge (n, rc) € E. U u,w € G, for some cluster G G Q, then ipc{u) 7 ^ ipc{w), and so 
7 ^ ^^{w). Otherwise, let Gu (respectively, Gw) be the cluster that contains u (resp., w), and Gu 7 ^ Gw 
The clusters Gu and Gw are adjacent in Q{Q), and thus label{Gu) 7 ^ label{Gw)- Hence label{u) 7 ^ label{w), 
and so V’(^) 7 ^ V'(rc). 

Note also that x(G(C')) < x(G), for every vertex subset G CV. The coloring if employs maxc'gQ{x(G(C'))}- 
. logn colors, i.e., 0(x(G') • • logn). □ 

We proved the following theorem: 

Theorem 6.2. For any n-vertex graph G = {V, E) and an integer parameter k = 1,2,..., an 0(r?l^ logre)- 
approximation of the optimal coloring for G can be computed in 0(3^ • /c^) time. 

In particular, by setting the parameter k to be an arbitrarily large constant we can get a distributed 
0 (n'^)-approximation algorithm for the coloring problem with a constant running time, for an arbitrarily 
small constant e > 0. (The running time is 0(31^^/*^^ • ^).) This greatly improves the current state- 
of-the-art constant-time distributed approximation algorithm for the coloring problem due to [5], which 
provides an approximation guarantee of 0(re^/^"*“^). On the other hand, the dependence of the running 
time on e is only 0(l/e) in the result of [5]. 

Note that the algorithm in Theorem 16.21 requires very heavy (exponential in n) local computations and 
large messages. The heavy computations are inevitable, because unless NP = P, the coloring problem 
cannot be approximated up to a ratio of for any constant e > 0 [271 [H Si- 

6.2 Coloring Triangle-Free Graphs and Graphs with Large Girth 

A result of Ajtai et al [ 2 ] shows that triangle-free n-vertex graphs G admit an 0{^/n|^/\ogn)-co\oxm.g. 

(This existential bound was shown to be tight by Kim [29]. i Here we show that one can construct an 
Q(j.ji/2-i-e)-coioring of triangle-free graphs in distributed randomized 0(l/e) time. Moreover, unlike our 
algorithm from the previous section, this algorithm uses only short messages and does not rely on heavy 
local computations. 

The algorithm starts with invoking the algorithm from Corollary 14.41 on its input n-vertex graph 
G = {V,E) with the parameter k = 2. We obtain a (2, 0(n^/^ log n))-network decomposition Q in 0(1) 
time. Moreover, the algorithm also constructs an Ff-partition ( 81 , 82 ) of the vertex set F of G into two 
sets. The degree of this FF-partition is A = logn). The clusters in 81 are singleton clusters. (Each 

such a cluster G G 81 contains a single vertex v G G such that deg(u) < A.) Each cluster G G ^2 is a star 
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rooted at a center vertex r G C. Also, since the graph is triangle-free, neighbors of r are not connected 
via edges one with another. 

Centers of clusters of S 2 now toss a color for their cluster from [A ■ rf]. If a color tossed by the root r of 
C is different from the colors of clusters incident on C in the supergraph Q{Q), then r stops. Otherwise it 
continues. Overall, as we have seen, after 0(l/e) rounds, clusters of S 2 will be colored in logn) 

colors. (The communication between centers of adjacent clusters can be executed efficiently using short 
messages. This requires some care. The root r of each cluster informs all vertices of C of its choice of 
color. Then each vertex of C (including r) sends the root’s color c(r) over inter-cluster edges incident on 
them. Then every vertex v & C checks if one of its neighboring clusters chose a color c equal to c(r). If 
it is the case, then it informs r. In this case r abandons its color (and informs all vertices of C about it), 
and continues to the next round of the randomized coloring procedure.) Then clusters of toss colors 
for them from the same range. Since each cluster of Si has only 0{A) neighbors in S'! U 5*2, the coloring 
will be computed within additional 0(l/e) rounds. Finally within each cluster C G S 2 actually two colors 
are used. (One for the center, and another for its neighbors.) Hence the overall number of colors is at 
most 2 ■ A ■ = 0{n^^‘^^^logn). The factor logn can be swallowed by slightly increasing the e in the 

exponent. To summarize: 

Theorem 6.3. An coloring of triangle-free n-vertex graph can he computed in 0(l/e) dis¬ 

tributed randomized time, using short messages and polynomially-bounded local computations. 

This result extends also to graphs with large girth. Specifically, consider a graph G with girth greater 
than g, for some integer g = 2k, k > 2. The arboricity of G is at most (See, e.g., [13] . Theorem 

3.7.) By [6], an H-partition S*!, S' 2 ,..., S'^, 1 = 0(l/e), of G with degree A = can be computed 

in constant time, for an arbitrarily small e > 0. Hence, by Theorem 14.51 an A ■ n^-coloring of G can be 
computed in 0(£/e) = 0(1) time. By scaling e we obtain the following result. 

Theorem 6.4. For a graph G with girth greater than g = 2k, k > 2, and an arbitrarily small constant 
e > 0, an coloring can he computed in constant distributed randomized time (specifically, 0(l/e^)j, 

using short messages and polynomially-bounded local computations. 

Note that the algorithm from Theorem 16.41 does not employ a network decomposition. Observe also 
that for k = 2 (i.e., girth greater than 4) the numbers of colors in Theorems 16.31 and 16.41 are the same, 
and both are existentially tight up to a slack factor of On the other hand, their proofs are different. 
However, Theorem 16.31 applies for g > 3 too, while Theorem 16.41 applies only for g > 4. So the result of 
Theorem 16.41 is mainly of interest for k>3 (i.e., g >()). 

6.3 Separated Decompositions 

For the sake of some applications we need a stronger notion of network decompositions, called a a- 
separated {a, l3)-network-decomposition, for positive parameters a, a, and fi [T]. An (a,/3)-network- 
decomposition Q of a graph G = (V, E) is called a-separated if the clusters of Q can be /3-colored in such 
a way that every pair of clusters G,G' ^ Q which are colored by the same color are at distance at least a 
from one another, i.e., distG(C', C*') > cr. Observe that an ordinary network decomposition is 2-separated. 

It is very easy to convert any procedure that constructs an ordinary (2-separated) (a, /3)-network- 
decomposition into a procedure that constructs a weak u-separated {a - (ct — 1 ), /3)-network-decomposition, 
for any parameter it > 3. (See Section [2] for the dehnition of weak decomposition.) Specihcally, one just 
executes the procedure for computing an ordinary {a, /3)-network-decomposition on the graph G^~^ = 
(y, = {(tt,u) \ u,v & V, distG(tf,u) < cr — 1}. As a result one obtains a partition Q of 

G^~^ such that each cluster G £ Q has diameter at most a in G^~^, and thus weak diameter at most 
(fj — 1) • a in G. Also, for any pair G, G' of distinct clusters in Q which are colored by the same color, the 
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distance between them in is at least 2, and so the distance between them in G is at least a. Hence 

Q is a weak ir-separated (a ■ (a — 1),/3)-network-decomposition of G. Simulating a distributed algorithm 
for G'^~^ in G increases the running time by a factor of a — 1. (Here we assume that message size is 
unbounded.) Therefore, Corollary 14.41 implies the following result. 

Corollary 6.5. For a pair of positive integer parameters cr > 2,k > 2, a a-separated weak ((3^“^ — 

1) • a, logn))-network-decomposition Q and an H partition {Si, S 2 , ■■■, Sk) of length k and degree 

log n) forQ{Q) can be computed in randomized 0{3^-a) time, with high probability. Moreover, for 
an integer parameter t = 1,2,one can compute an 0{r?l^ log^ nlog^*^ n)-labeling for Q in 0{3^ ■ a - t) 
time. 

We remark that this simple approach for converting network-decompositions into weak separated ones 
is not new. It was used, e.g., by Dubhashi et al. [El- 

Next we show that our algorithm for constructing ordinary (3^ ^ — 1, 0{n^/^ log n))-network-decompositions 
can be adapted to compute strong cj-separated ((2cr — — 1, 0{n^^^ logn))-network-decomposition in 

randomized time 0{{2a)^), for an arbitrary integer parameter cr > 2. 

In what follows we describe Procedure Sep-Decompose which generalizes Procedure Decompose (Al¬ 
gorithm [1]). It accepts as input all the parameters of Procedure Decompose, and also the separation 
parameter a. Consider again Procedure Decompose (Algorithm [T]) . The termination condition of the 
procedure (lines 1-2, the case when the size s is small, i.e., s = 0(n^/^ log n)) stays unchanged. In the 
general case (the ’’else” case of the procedure, lines 3-15) Procedure Decompose starts with invoking 
Procedure Partition, which decomposes the input graph G into A and B. In the original procedure the 
subgraph G{A) induced by A has a small maximum degree (at most 0{qlogn), where q = 0(n^/^) is 
an input parameter of Procedure Partition.) The generalized variant of the procedure (Procedure Sep- 
Decompose) invokes instead a generalized variant of Procedure Partition, called Procedure Sep-Partition. 

The latter procedure accepts as input all the parameters of Procedure Partition, but also the separation 
parameter a. It also decomposes the input graph G into A and B, but A has the property that 
has maximum degree 0{q\ogn) = 0{n^^^\ogn), i.e., for every vertex u E A, there are at most 0{q\ogn) 
other vertices of G at distance at most a — 1 from v. (The distance is with respect to G.) Similarly to 
Procedure Partition, in Procedure Sep-Partition too the set H is a collection of at most sjn^^^ clusters of 
small diameter in G. However, the diameter grows from 3 in the case of Procedure partition, to 2(T — 1 
in Procedure Sep-Partition. 

Then Procedure Sep-Decomposition invokes Procedure Dec-Small. (See line 5 of Algorithm [TJ) Pro¬ 
cedure Dec-Small converts every vertex C € A into a separate cluster. (If G is the original graph G then 
every vertex u E A C P gives rise to a cluster {u}. Otherwise G is a supergraph of the original graph G, 
and a vertex C E A is a cluster of G.) The resulting set of clusters is denoted by S. Procedure Dec-Small 
also returns a labeling for clusters of S, but similarly to the case of Section 01 this labeling is immaterial 
for our discussion. 

On line 6 of algorithm [1] Procedure Sep-Decompose invokes itself recursively on the supergraph Q{B) 
induced by the set B of clusters. The rest of the procedure stays unchanged. 

At this point we are interested in a version of Procedure Sep-Decompose which only computes a 
separated network-decomposition without a labeling function for it. (See the beginning of Section 01) 

To recap, this procedure returns a network-decomposition Q = where {Si, S 2 ,..., Sk) is an H- 

partition of the supergraph Q{Q) induced by this decomposition. (See Lemmas 0A] and 021) Moreover, it 
is easy to verify that decompositions Q produced by Procedure Sep-Decompose satisfy a stronger property 
than decompositions produced by Procedure Decompose. Specifically, by construction, for every index 
i = 1,2, ...,k, a cluster C ^ Si has at most 0{q\ogn) = 0(n^/^ logn) other clusters G' E at 

distance at most cr — 1 from it in G. This fact is summarized in the next lemma. 
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Lemma 6.6. {Si, S 2 , ■■■, S^) is an H-partition with degree 0(n^/^logn) of the supergraph {G{Q))^ 
where Q = U^^iSi. 

By invoking one of the algorithms from Section |3] for coloring low-arboricity graphs (for which a short 
low-degree //-partition is provided) we can obtain an log^ nlog^*^ n)-labeling for Q, which has the 

property that any two distinct clusters C, C which receive the same label are at distance at least a from 
one another in Q{Qf), and thus at distance at least a from one another in G. The running time of this 
step is 0{t ■ Diam{Q)). (See Theorem 14.41 1 Alternatively, one can have an 0(n^/^ log n)-labeling with 
this property in time 0{Diam{Q) ■ k'^). (See Corollary 14.71 1 

Next we analyze Diam{Q). To do it we hrst describe Procedure Sep-Partition. (See Algorithm [3] for 
Procedure Partition.) The procedure accepts the same parameters as Procedure Partition, but also the 
separation parameter a. (In fact. Procedure Partition is a special case of Procedure Sep-Partition, where 
a = 2.) Similarly to Procedure Partition, in Procedure Sep-Partition every vertex v selects itself (joins 
D) independently at random with probability 1/q. Then every selected vertex v sends an exploration 
message to distance cr — 1 from it in G. Every vertex u which is not selected {u ^ D) and receives 
at least one exploration message joins the cluster centered by the closest originator of an exploration 
message received by u. (Ties are broken in an arbitrary but consistent way according to the identities of 
originators. If originators themselves are clusters, then each of them has its own leader whose identity 
serves as the identity of the cluster. The consistent rule for breaking ties may be, for example, to prefer 
an originator with a smaller identity.) Other vertices join the set A. The procedure returns the set A 
and the set B of clusters which are created in the way described above. Observe that if G is not the 
original graph but rather a supergraph of it then the algorithm is executed by clusters rather than by 
single vertices. In other words, in this case the center of each cluster simulates all the operations that 
need to be performed by the cluster. 

The next lemma shows that clusters created by Procedure Sep-Partition are connected and have 
bounded diameter. 

Lemma 6.7. Consider an invocation of Procedure Sep-Partition(G,q,a), where q > 1 is a parameter 
and a >2 is an integer parameter. Then each vertex v G A has degree 0{q ■ logn) in and each 

cluster G & B has (strong) diameter at most 2a — 2 in G. 

Proof. Let c be a sufficiently large fixed constant, and consider a vertex u G G such that a {a — 1)- 
neighborhood Bg— i(u) of u in G contains at least c ■ q ■ logn vertices. Then with probability at least 
1 — 1/n^ at least one of the vertices u G i?(j-i(u) joins D, and the vertex v becomes clustered. Hence 
with probability at least 1 — l/n^“^ all vertices v with |Ho-_i(u)| > c ■ q ■ logn become clustered, and so 
each unclustered vertex v G A satisfies \B„-i{v)\ < c - q ■ logn. 

Consider a cluster G G B. It is centered around an originator v of an exploration message. (The 
vertex v belongs to D, i.e., it is selected.) Consider a vertex u G G, and let P^^u be a shortest v — u path 
in G. Let x be a vertex on this path. (Note that v,u,x are vertices of G, i.e., they are possibly clusters 
themselves.) It follows that v is the closest selected vertex to x, and if there exists another selected vertex 
v' G D which satisfies distg,(u,x) = distg,(u', x), then v has a smaller identity than v'. (As otherwise v' 
would rule u as well.) Hence x G C. Consequently all vertices of are in G, and the length of Pu^v 
is at most a — 1. Hence the cluster G has strong radius at most a — 1, i.e., strong diameter at most 
2{a-l). □ 

Observe also that by the same argument as in Lemma 13.101 the number of clusters in B is, with 
high probability, 0{s/q). We are now ready to analyze the diameter Diam{Q) of the ultimate network- 
decomposition Q. The following lemma generalizes Lemma 13.31 
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Lemma 6.8. Let Q = (^i, 5*2 ,5^) he a a-separated network-deeomposition produced hy the invocation 
Sep-Decompose(G,n,k, s := n,e,t,a) on an input graph G. Then for each i E [k], Diam{Si) < {2a — 
l)*-i - 1. 

Proof. We prove by induction on i that in the rth level recursive invocation of Procedure Sep-Decompose 
each vertex v of the input graph Gi of this invocation is a cluster of the original graph G with diameter 
at most {2a — 1)*“^ — 1. Since for each i E [A;], clusters of Si are vertices of Gi = {Vi,Ei) the assertion of 
the lemma follows from the inductive claim. 

Base: Diam{Si) = 0 = {2a — 1)*^ — 1. 

Step: Consider an index i < k. By Lemma 16.71 each cluster G created by the ith level invocation of 
Procedure Sep-Decompose has strong diameter at most 2a — 2 in Gi = {Vi,Ei). It follows that 

Diam{G) < {2a — 1) • max{Diam(C")} -|- {2a — 2). 

C'eVi 

By induction hypothesis it follows that 

Diam{G) < {2a - l){{2a - 1)*"^ - 1) + {2a - 2) = {2a - If - 1. 

Since vertices v of Gj+i are clusters which were formed by the ith level invocation of Procedure Sep- 
Decompose, the assertion of the lemma follows. □ 

We summarize this discussion with the following corollary. 

Corollary 6.9. Consider an invocation of Sep-Decompose(G,n,k,s := n,e,t,a), where k > l,a >2 are 
integer parameters. It produces a a-separated strong {{2a—l)^~^ — l,0{n^^^logn))-network-decomposition 
Q = U^^■^^Si, along with an H-partition (S'i,52, ...,Sk) for {Q{Q)y~^. The running time of this invocation 
is 0{{2a — 1)^“^). 

As was discussed in the paragraph following Lemma 16.61 using this network-decomposition one can 
compute an 0{nfl^ log^ nlog*-*^ n)-labeling for Q within additional Oft ■ {2a — rounds, or alterna¬ 

tively, an Ofn?/^ log n)-labeling within additional 0{{2a — 1)^~^ ■ k'^) rounds. In both cases the labeling 
satisfies that any two distinct clusters G, G' which receive the same label are at distance at least a one 
from another in G. 

One can also improve the parameters of the network-decomposition from Corollary 16.91 from ((2 ct — 
1)^“^ — 1, Ofn}/^ log n)) to {0{a)^,0{n^^^)) at the expense of increasing the running time from {{0{2a — 
1)^“^) to log^/^ n in general graphs, and 0{a)^ ■ exp{0{y/\og log n)} in graphs with girth at least 

6. This is done by introducing to Procedure Sep-Decompose a modification analogous to the one that 
we introduced to Procedure Decompose in Section [5l Recall that the difference between Procedure RS- 
Decompose and Procedure Decompose is that the former invokes Procedure RS-Partition as a subroutine, 
while the latter invokes Procedure Partition. 

Procedure RS-Partition computes a (3,5)-ruling set W for the set [/ = {ri E P | deg(n) > q} of high 
degree vertices of its input graph G, for a parameter 6. The variant of this procedure that we are now 
describing, called Procedure Sep-RS-Partition, accepts as input also the separation parameter a, and 
computes a {2a — l,5)-ruling set W for the set U' = {u ^ U : |i?o-_i(n)| > q} of vertices that have 
at least q vertices in their {a — l)-ball. The clusters {Gw \ w E W'{ are then created in the same way 
as in Procedure RS-Partition. In particular, their strong radii are still bounded by 5. Also, every vertex 
u ^ U' is assigned to some cluster. The sets A and B are now formed as in Procedure RS-Partition. 
Every vertex u E A now satisfies \B„-i{v)\ < q. The following lemma is analogous to Lemma l5.11 and its 
proof is very similar to that of Lemma 15.11 
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Lemma 6.10. Suppose that Procedure Sep-RS-Partition is invoked on a graph G = {V,E) and positive 
parameters a and q. Suppose further that it uses a subroutine for computing a {2a — l,6)-ruling set, for 
a positive integer parameter a. Then in the subgraph G{A) every vertex v ^ A satisfies \B^_i{v)\ < q. 
Moreover, B consists of at most \ V\/q supernodes, each of which is a cluster of strong diameter at most 
25. 

It follows now that Procedure Sep-RS-Decompose computes a cr-separated ((2(5+1)^, 0(n^/^))-network- 
decomposition Q. For the running time we need again to specify the running time required for computing 
a (2cr — l,(I)-ruling set. By running the algorithms for computing a ruling set due to Barenboim et al. 
m and Kothapalli and Pemmaraju m respectively in G'^^^ we obtain a {2a — 1 , 2 - 2 cj)-ruling set in 
the case of general graphs, and a {2a — 1,3- 2cr)-ruling set in the case of graphs of girth at least 6 . In both 
cases 5 = 0{a), and the running time is 0{a ■ log^^^n) in the former case and 0{a) ■ exp{0{^/\og\ogn)} 
time in the latter. The rest of the analysis is identical, except that the overall running time of Procedure 
Sep-RS-Decompose becomes (0 (cj))^ log^^^ n and {0{a))^ ■ exp{0{ \/log log n)} in the cases of general 
graphs and graphs of girth at least 6 , respectively. 

Theorem 6.11. Procedure Sep-RS-Decompose invoked on an n-vertex graph G = {y,E) with positive 
integer parameters k and a computes a a-separated strong {{0{a))^,n^^^)-network-decomposition Q in 
randomized time {0{a))^ n in general graphs and in {0{a))^ ■ exp{0{^J\og\ogn)} randomized time 

in graphs of girth at least 6 . 

One application of strong separated network-decomposition is low-intersecting partitions. Low-intersecting 
partitions were introduced by Busch et al. in their work on universal Steiner trees. A low-intersecting 
{a, -partition V of a graph G is the partition of the vertex set V such that 

(1) Every cluster G in V has strong diameter at most a • 7 . 

(2) For every vertex u E R, a ball B-^{v) of radius 7 around v intersects at most /3 clusters of V. 

Busch et al. showed that given a hierarchy of low-intersecting partitions with certain properties (see 
m for details) one can construct a universal Steiner tree. (See [13] for the definition of universal Steiner 
tree.) Also, vice versa, given universal Steiner tree they showed that one can construct a low-intersecting 
partition. They constructed a low-intersecting partition with a = 4^,/3 = k ■ and arbitrary 7 . 

We next argue that a ( 27 -|-l)-separated strong {p, 77 )-network-decomposition Q is also a low-intersecting 
partition with parameters {a = /i/ 7 ,/3 = 7 , 7 ). Indeed, every cluster G of Q has strong diameter at most 
/i = a • 7 . Moreover, consider a vertex v and a ball B-y{v) of radius 7 around v. Observe that for every 
color class i E [r/] of G{Q), the ball B.y{v) can intersect at most one cluster G colored by i. (This is 
because for every two i-colored clusters G,G', it holds that distG(C', C") > 27 - 1 - 1 .) Hence altogether 
B-f{v) may intersect up to r/ clusters of Q. This proves the claim. 

Therefore, our distributed algorithm for computing a ( 27 -|-l)-separated strong ( 0 ( 7 ))^, n^/*^)-network- 
decomposition in distributed randomized time ( 0 ( 7 ))^ log^^^ n in general graphs and in 
( 0 ( 7 ))^ • exp{0{\/log logn)} in graphs of girth at least 6 provides also a distributed algorithm with the 
same running time for constructing a low-intersecting (( 0 ( 7 ))^, 7 )-partition. We summarize: 

Corollary 6.12. For any pair of positive integer parameters A:, 7 , a low-intersecting (( 0 ( 7 ))^, 7 )- 

partition can be constructed in ( 0 ( 7 ))*^ log^^^ n randomized time in general graphs and in ( 0 ( 7 ))^ • 
exp{0{y/\og\ogn)} randomized time in graphs of girth at least 6. 

We remark that this construction can be implemented using short messages. 

Comparing this result with the algorithm of Busch et al. m we note that the partition of m has smaller 
radius. (It is 7 - (0(1))^ instead of ( 0 ( 7 ))^ in our case.) On the other hand, the intersection parameter /3 
of our partitions is smaller. (It is instead of k-n^^^.) In particular, the intersection parameter in the 
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construction of [13] is always n(logn), while ours can be as small as one wishes. Finally, the algorithm 
of Busch et al. m is not distributed, and seems inherently sequential. 


6.4 Approximation Algorithms for the Minimum Dominating Set and Minimum t- 
Spanner Problems 


In this section we employ our network-decomposition algorithm in order to derive approximation algo¬ 
rithms for the minimum dominating set and minimum t-spanner problems. Suppose that we are given a 
3-separated (d, £)-network-decomposition Q of a graph G. For each cluster C £ Q, we compute in paral¬ 
lel a dominating set D C F'’'(C) of C, such that D has minimum cardinality among all dominating sets 
D' C F'’'(C) of C. The computation of D is performed by collecting the topology of the clusters and their 
neighborhoods by the leaders of respective clusters, performing the computation locally using exhaustive 
searclJi] , and broadcasting the results to the vertices of the clusters and their neighbors. Since the weak 
diameter of the clusters is at most d, this requires 0{d) rounds. We next show that the resulting set 
obtained by taking the union of the dominating sets in all clusters constitutes an .^-approximate minimum 
dominating set of the input graph G. 

Lemma 6.13. For a 3-separated {d, i)-network-decomposition Q, suppose that we have computed a mini¬ 
mum dominating set Dq F r'''(C) of C, for each cluster G G Q. Then \ \ C G Q}\ < £-\MDS{G)\. 

Proof. For 1 < i < i, let Ui O V denote the set of all vertices with label i in the network-decomposition Q. 
Let lli = r'*'(C/j). We claim that | IJi^c I C C Ui}\ < \W\, where VF is a minimum dominating set of G. 
(Note that in the current proof the notation G ^Ui stands for a cluster G that belongs to C/j, rather then 
just a subgraph of Ui, since Dq is defined only for clusters.) Let Gi £ Q he a, cluster of label i, 1 < i < £. 
Then C* C Ui. Observe that W n r+(Cj) is a dominating set of Cj. (Since IF is a dominating set of Cj, 
and any vertex in IF \ F'’'(Ci) does not dominate any vertex in Cj.) Therefore, |IF n r'^(C'i)| > \Dcf\. 
Note also that for any cluster C[ / Ci of label i it holds that r+(C'') n r+(C'j) = 0. Indeed, Q is a 
3-separated network-decomposition, and thus, for any u £ Ci,v £ C[ it holds that distG(n, n) > 3. Hence 
for any w £ F+(C'j),x £ r+(C''), it holds that distG'(a:, > 1, and thus x w. Consequently, 


\^{Dc\ GCUi}\ = E i»ci < E |IFnF+(Ci)| < |1F| = \MDS{G)\. 

CCUi CQUi 

Therefore, 

e 

I |J{Zlc I c G Q}| = I U {{U{Dc \G<ZUi})-.i£ [£]} I < I U {He I c C C/JI < £ • \MDS{G)\. 

i=\ 

□ 


Recall that by Corollary 15.41 there is a routine that computes an ((0(logn^/^)-network- 
decomposition in deterministic time (0(logn))^“^, for any k = 1 , 2 ,.... As was discussed above, this 
routine can also be adapted to compute a weak 3-separated network-decomposition with the same prop¬ 
erties within the same running time. (See Section [6.31 both the diameter parameter and the running time 
grow by a constant factor a = 3.) Also, similarly to Corollary 16.91 one can adapt this routine so that 
it will compute a strong network-decomposition with the same parameters and the same running time. 
(The diameter and the running time grow by a factor of (2cj — 1)^ = 5^, which is however swallowed 
by the notation (0(logn))^“^.) Using this network-decomposition in conjunction with Lemma 16.131 we 
obtain the following theorem. 

^We note that once can employ polynomial-time local computations instead of exhaustive search in the expense of 
increasing the approximation ratio by a factor of 0(log A). See Section [T] 
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Theorem 6.14. For an n-vertex graph G, and a positive integer parameter k an approximation 

for the minimum dominating set problem can he computed in deterministic time ( 0 (logn))^“^. 

Another problem for which an efficient approximation algorithm can be obtained using network- 
decompositions is the minimum t-spanner problem. Given an (unweighted) graph G = (1/, E) and a 
positive integer parameter t, a subgraph G' = H C E, is a t-spanner of G if for every pair 

u,v G V of vertices, distc'(u,u) < t ■ distcCw,t')- In the minimum t-spanner problem the objective is to 
find a t-spanner of the input graph with as few edges as possible. 

Suppose that we are given a ( 2 t — l)-separated (d, .^)-network-decomposition Q of an input graph 
G = {V, E). Let Cl, C 2 ,..., Gh be a single color class of this labeling, i.e., Diam{Gi) < d for each i G [/i], 
and distG(Ci, Gj) > 2t — 1, for every pair of distinct indices i 7 ^ j, i,j G [h]. Let Cj = Bt-i{Gi), for every 
i G [h]. Note that Cj n Cj = 0, for every pair of distinct indices i / j. Denote C = u(LiCj,C = U^LiCj, 
and consider a minimum t-spanner Ti for E{C) which is allowed to use edges from E{C). Let also H* be 
a minimum t-spanner for G. 

Lemma 6.15. \FL\ < \H*\. 

Proof. Observe that the restriction H*{C) of H* to C is a t-spanner for C. Indeed, consider an edge 
{u,v) G E{C). Let Gi G C he the cluster such that u,v G Ci. Then H* contains a path of length at most 
t between u and w, and so this path belongs to H*{Ci) C H*{C). 

The lemma now follows as T-L is the minimum t-spanner for C which is allowed to use edges from E{G), 
while H*{C) is a t-spanner for C of this type. Hence \'H\ < \E[*{C)\ < \E[*\. □ 

Denote also by Hi the minimum t-spanner for E{Gi) which is allowed to use edges of E{Gi). 
Lemma 6.16. | U^^i Hi\ = \T-L\. 

Proof. Obviously, is a t-spanner for E{C) which uses only edges of E{C). Hence by optimality of 

In the opposite direction, for every index i G [h], let Hi = HFE{Ci). By optimality of Hi, \Hi\ < \Hi\. 
Also, for every pair of distinct indices i,j G [h], HifiHj = 0 . (This is because E{Gi)d\E{Gj) = 0 .) Hence 

h h 

\n\ = I uti H,\ = 1^*1 > E 1^*1 ^ I uti H^\. 

i=l i=l 

(The last inequation is, in fact, equality.) □ 

In other words, to compute a minimum t-spanner H for E{C) one can compute minimum t-spanners 
Hi, H 2 ,..., Hi for E{Gi), E{C 2 ), ..., E{Gh) (which are allowed to use edges of E{Ci), E{C 2 ), ■■■, E{Ck), 
respectively), and take their union. Our distributed algorithm will do precisely this. In each cluster G 
of Q it computes a minimum t-spanner for E{C) using edges of E{C), C = Bt{G). This computation is 
done by collecting the entire topology of {C, E{G)) into a vertex in G, doing a local (possibly very heavy) 
computation, and informing all vertices of G about the results of this computation. The union of all these 
t-spanners will be our ultimate spanner. Hence the algorithm returns a spanner H' = where 

for each index j G [h\, H^^^ is a minimum t-spanner for E{C^^'>), where is the set of all vertices labeled 
by j in the network decomposition Q. (In other, words, they belong to clusters of color j. Note, however, 
that to execute the algorithm we do not need to know these colors/labels.) Since by Lemma 16.161 for 
every j G [i\, < \H*\, it follows that the algorithm returns an ^-approximation. The running time 

of the algorithm is 0{Diam{Q) + t) = exp{0{k)} + 0{t). To summarize: 
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Theorem 6.17. For any pair of positive integer parameters t,k an 0{n^^^)-approximation of the mini¬ 
mum t-spanner problem in n-vertex graphs can be computed in exp{0{k)} + 0{t) randomized time. 

Observe that the same result applies to the t-spanner problem in directed graphs, by the same argu¬ 
ment. Note that even though the graph is directed, we assume that the communication over every edge 
is bidirectional. 

7 Removing heavy local computations from the minimum dominating 
set and minimum t-spanner algorithms 

It is well known that an 0(log A)-approximation of minimum dominating set can be computed in poly¬ 
nomial time in the sequential setting. (See, e.g., [50].) However, this approach cannot be applied directly 
to our algorithms since we compute minimum dominating sets Dq of clusters C, such that Dc C r'’'(C') 
rather than Dq C C. On the other hand, this problem reduces to the Set Cover problem with both the 
degree parameters (i.e., the maximum cardinality of a set and the maximum number of sets that share an 
element) bounded by A -|-1. Hence this problem admits a polynomial-time 0(log A)-approximation algo¬ 
rithm. (See, e.g., [l 6 |.) One can also extend the classical centralized 0(log A)-approximation algorithm 
for the MDS problem directly to our slightly more general problem. This extension is described below. 
Consequently, we can obtain a dominating set whose size is at most O(logA) the size of the minimum 
dominating set of C consisting of vertices of r'*'(C'). This can be achieved in the following way. Initially 
Dc = 0. We proceed in phases, each time selecting a vertex v from C such that dy = |r'''(u)nC'\r’^(T)c)| 
is maximal, and adding v to Dc- (Ties are broken by preferring vertices that belong to T~^(Dc), and if 
this does not solve the tie, it is broken arbitrary.) Once no uncovered vertex remains we are done. 

Let S* C r'*'(C') be a minimum dominating set of C. We claim that Dc < 0(log A • l^*!). The proof 
is by amortized analysis. Each time a vertex v is added to Dc we assign a weight 1/d^ to each vertex of 
r+(u) n C \ r+(L>c'). Observe that the sum of all weights assigned during this procedure is \Dc\. Next, 
(for analysis) let each vertex of C select a single vertex from S* that dominates it. Consider a vertex 
u G S* and the set W of all neighbors of u in C that selected u. Next, we analyze the sum of weights of 
W. For each tc G IT it is assigned a weight once a neighbor of w (or w itself) joins the dominating set 
Dc- Let i be the number of the phase in which it happens, and degj(ri) denote the number of neighbors 
of M in C that are not covered in the beginning of phase i. Also, let z denote the neighbor that dominates 
w, for which w obtained its weight. Since in each phase a vertex v with maximal is selected, it holds 
that degj(ri) < degj(z). Consequently, w is assigned a weight at most 1/ degj(z) <1/degj(ri). Therefore, 
the sum of weights of W is at most 1/j = 0(log A). Therefore, the sum of all weights in the graph 

is iZ^cl = 0(log A • I S'*!). This completes the proof. 

A similar idea can be applied in the case of the minimum t-spanner problem. Again, we need a 
centralized polynomial-time approximation algorithm for the minimum t-spanner for edges of E{C) (for 
a cluster (7), while the spanner is allowed to use edges of E{C). This is an instance for the client-server t- 
spanner problem, and for the case t = 2 it was devised in |2T] . By plugging it in our distributed algorithm 
for approximating spanners we obtain a distributed 0 (n^/^ log n)-approximation algorithm with running 
time exp{0{k)'\ for the directed and undirected 2-spanner problem. The latter algorithm only employs 
polynomially-bounded local computations. To the best of our knowledge, there are no existing centralized 
algorithms with a non-trivial approximation guarantee for the client-server t-spanner problem for t > 3. 
It is however likely that the LP-based approaches to the minimum t-spanner problem (such as [111 [T8] ) 
extend to this more general problem. 
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